#include <gmlib/number/mont256.h>
#include <gmlib/number/uint256.h>

#include <stdexcept>
#include <cstring>

using namespace number;

static std::uint8_t  P_data[32]      = {0xff,0xff,0xff,0xfe,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff};
static std::uint8_t  P_SUB2_data[32] = {0xff,0xff,0xff,0xfe,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0xff,0xff,0xff,0xfd};
static std::uint8_t  R_data[32]      = {0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01};
static std::uint8_t  R_POW2_data[32] = {0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x02,0xff,0xff,0xff,0xff,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03};
static std::uint8_t  N_data[32]      = {0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01};

void test_mont256_mul()
{
    uint256_t P, P_SUB2, R, R_POW2, N_;
    uint256_from_bytes(P, P_data);
    uint256_from_bytes(P_SUB2, P_SUB2_data);
    uint256_from_bytes(R, R_data);
    uint256_from_bytes(R_POW2, R_POW2_data);
    uint256_from_bytes(N_, N_data);

    Mont256CTX ctx;
    ctx.P = P, ctx.P_SUB2 = P_SUB2;
    ctx.R = R, ctx.R_POW2 = R_POW2;
    ctx.N_ = N_;

    mont256_t    multiplier, multiplicand, product;
    std::uint8_t product_data[32];

    std::uint8_t multiplier0[32]   = {0xd2,0x3f,0x08,0x24,0x12,0x8b,0x2f,0x33,0x0c,0x5c,0x7f,0xd0,0xa6,0xa3,0xa4,0x50,0x65,0x13,0x27,0x0e,0x26,0x9e,0x0d,0x37,0xf2,0xa7,0x4d,0xe4,0x52,0xe6,0xb4,0x38};
    std::uint8_t multiplicand0[32] = {0x36,0xf6,0x75,0xcc,0x81,0xe7,0x4e,0xf5,0xe8,0xe2,0x5d,0x94,0x0e,0xd9,0x04,0x75,0x95,0x31,0x98,0x5d,0x5d,0x9d,0xc9,0xf8,0x18,0x18,0xe8,0x11,0x89,0x2f,0x90,0x2b};
    std::uint8_t product0[32]      = {0xf2,0x8e,0x23,0x8a,0xcf,0x20,0x77,0x3b,0xa1,0x36,0xe6,0xfc,0xdd,0xc4,0x78,0xa2,0x21,0x1c,0x89,0x39,0xc4,0xbd,0x40,0xa6,0x91,0x6e,0x48,0x04,0x52,0xf3,0x3e,0x30};
    mont256_from_bytes(&ctx, multiplier, multiplier0);
    mont256_from_bytes(&ctx, multiplicand, multiplicand0);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product0, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier1[32]   = {0x8d,0x11,0x6e,0xce,0x17,0x38,0xf7,0xd9,0x3d,0x9c,0x17,0x24,0x11,0xe2,0x0b,0x8f,0x6b,0x0d,0x54,0x9b,0x6f,0x03,0x67,0x5a,0x16,0x00,0xa3,0x5a,0x09,0x99,0x50,0xd8};
    std::uint8_t multiplicand1[32] = {0xa1,0x70,0xb3,0x38,0x39,0x26,0x30,0x59,0xf2,0x8c,0x10,0x5d,0x1f,0xb1,0x7c,0x23,0x90,0xc1,0x92,0xcf,0xd3,0xac,0x94,0xaf,0x0f,0x21,0xdd,0xb6,0x6c,0xad,0x4a,0x26};
    std::uint8_t product1[32]      = {0x75,0x0f,0x02,0xa1,0x4e,0x2a,0xd3,0xf7,0xbf,0x5f,0x0d,0xc2,0xc0,0x7b,0xd5,0x32,0x01,0xde,0x10,0x68,0x5f,0x3b,0x0f,0xbf,0x02,0x6a,0x50,0x49,0xaa,0x4b,0x01,0x9d};
    mont256_from_bytes(&ctx, multiplier, multiplier1);
    mont256_from_bytes(&ctx, multiplicand, multiplicand1);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product1, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier2[32]   = {0x0c,0xb1,0xe2,0x9c,0x65,0x8c,0xda,0x14,0x95,0xe6,0x0a,0xf5,0x93,0xbd,0x04,0xcf,0x0f,0xd6,0x30,0xf1,0xf2,0x9d,0x0d,0xa9,0x95,0x3f,0x48,0xf1,0xa0,0x9f,0x76,0xb5};
    std::uint8_t multiplicand2[32] = {0x6b,0x4c,0xb2,0x42,0x4a,0x23,0xd5,0x96,0x22,0x17,0xbe,0xad,0xdb,0xc4,0x96,0xcb,0x8e,0x81,0x97,0x3e,0x0b,0xec,0xd7,0xb0,0x38,0x98,0xd1,0x90,0xf9,0xeb,0xda,0xcc};
    std::uint8_t product2[32]      = {0x2a,0x34,0x5a,0x50,0xd5,0xd7,0x3d,0x30,0x0c,0xf0,0xea,0x59,0x78,0xb9,0x75,0x0a,0x5a,0x25,0xa6,0x3e,0xb8,0xc7,0xfd,0x66,0xd1,0x60,0xb2,0xc4,0x09,0x31,0xc2,0x52};
    mont256_from_bytes(&ctx, multiplier, multiplier2);
    mont256_from_bytes(&ctx, multiplicand, multiplicand2);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product2, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier3[32]   = {0xae,0x97,0xba,0x94,0xd0,0xed,0xa8,0x2f,0x8f,0x6d,0x05,0x58,0x4e,0xf8,0xaa,0x38,0x92,0x27,0x66,0x58,0x1e,0x27,0xa1,0xc0,0x8a,0x6a,0x63,0xec,0x24,0xed,0xe6,0xa4};
    std::uint8_t multiplicand3[32] = {0x18,0xf1,0x35,0xd2,0x5f,0x55,0x72,0x03,0x30,0x18,0x50,0xc5,0xa3,0x8f,0xd5,0x47,0x92,0x3a,0x73,0x69,0x94,0xe3,0xbf,0x91,0x1a,0x61,0xdb,0xe2,0x2e,0x44,0x15,0x8b};
    std::uint8_t product3[32]      = {0x79,0x12,0x44,0xd4,0x6d,0xfa,0xc6,0x9a,0xf9,0x4f,0x84,0x90,0x58,0xfc,0x05,0xbc,0x57,0xee,0xa7,0x5f,0xa8,0x78,0x94,0x0d,0x82,0x9e,0xc0,0xa1,0x58,0x24,0x8d,0xd2};
    mont256_from_bytes(&ctx, multiplier, multiplier3);
    mont256_from_bytes(&ctx, multiplicand, multiplicand3);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product3, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier4[32]   = {0x7f,0x15,0x05,0x24,0x34,0xb9,0xb5,0xdf,0x9e,0x77,0x69,0xb1,0x0f,0x42,0x05,0xb4,0x90,0x7a,0x70,0xc3,0x10,0x12,0xf0,0x37,0xb6,0x4c,0xe4,0x22,0x8c,0x38,0xfb,0x29};
    std::uint8_t multiplicand4[32] = {0xec,0x66,0xa7,0x87,0x95,0xe7,0x61,0xd1,0x77,0x31,0xaf,0x10,0x50,0x6b,0xf2,0xef,0xc6,0xf8,0x77,0x18,0x6d,0x76,0xb0,0x7e,0x88,0x1e,0xd1,0x62,0xae,0x2e,0xb1,0x54};
    std::uint8_t product4[32]      = {0x0f,0x98,0x55,0x98,0x15,0x7f,0x6e,0xf0,0x48,0x54,0x26,0x10,0xcb,0xbb,0xeb,0x0e,0x15,0x77,0xa4,0x96,0x06,0x0f,0x78,0xbc,0x84,0x5c,0x21,0x89,0x9c,0x64,0x50,0x0d};
    mont256_from_bytes(&ctx, multiplier, multiplier4);
    mont256_from_bytes(&ctx, multiplicand, multiplicand4);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product4, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier5[32]   = {0xc7,0xa2,0xea,0x20,0xb2,0xf1,0x4c,0x94,0x2e,0x05,0x31,0x9a,0xcb,0x5c,0x74,0x27,0x3f,0x98,0xe2,0x77,0x4c,0xbd,0x87,0xad,0x5c,0x90,0xa9,0x58,0x74,0x03,0xe4,0x30};
    std::uint8_t multiplicand5[32] = {0x57,0xee,0x05,0xcd,0xe0,0x09,0x02,0xc7,0x7e,0xbf,0xf2,0x06,0x86,0x73,0x47,0x21,0x4c,0xdd,0x20,0x55,0x93,0x0d,0x6e,0xaf,0x14,0xf4,0x73,0x3f,0x3e,0x7d,0x1b,0xfb};
    std::uint8_t product5[32]      = {0xaf,0xb4,0x82,0x28,0x49,0x4a,0x10,0x2e,0xf3,0x02,0x79,0xb5,0x6e,0x8e,0xbc,0x12,0x12,0x29,0x62,0x84,0xcf,0xa6,0x3e,0x2d,0x4f,0xf3,0xf5,0xa6,0x51,0xe6,0x89,0xb4};
    mont256_from_bytes(&ctx, multiplier, multiplier5);
    mont256_from_bytes(&ctx, multiplicand, multiplicand5);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product5, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier6[32]   = {0x83,0x0e,0x07,0xbc,0x1e,0x39,0x8f,0x10,0x12,0xbd,0x4a,0xce,0xfa,0xec,0xbd,0x38,0x9b,0xe4,0xbc,0xfc,0x49,0xb6,0x4a,0x08,0x72,0xe6,0xcc,0x3a,0xba,0xbc,0xed,0x20};
    std::uint8_t multiplicand6[32] = {0x6b,0xf4,0x6c,0x69,0x7d,0x2c,0xaf,0x82,0xee,0xea,0xcb,0xe2,0x26,0xe8,0x75,0x55,0x57,0x90,0xf8,0x2e,0xc1,0xd3,0xfc,0xff,0x2a,0x3a,0xf4,0xd4,0x6b,0x0a,0x18,0xe8};
    std::uint8_t product6[32]      = {0xdf,0x38,0xca,0x15,0x47,0x39,0xc9,0xee,0x1c,0xc7,0x8d,0x45,0x89,0x35,0x4c,0xe3,0xe3,0x93,0xc0,0x72,0x28,0x5e,0x5c,0x83,0x3b,0x80,0x01,0xe5,0x60,0xe2,0xbc,0x3c};
    mont256_from_bytes(&ctx, multiplier, multiplier6);
    mont256_from_bytes(&ctx, multiplicand, multiplicand6);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product6, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier7[32]   = {0xca,0x02,0x13,0x5e,0x92,0xb1,0xd3,0xf2,0x8e,0xde,0x0d,0x7a,0xc3,0xba,0xea,0x9e,0x13,0xde,0xef,0x86,0xab,0x10,0x31,0xd0,0xf6,0x46,0xe1,0xf4,0x0a,0x09,0x7c,0x97};
    std::uint8_t multiplicand7[32] = {0x7f,0x26,0x14,0x4b,0x98,0x28,0x9f,0xcd,0x59,0xa5,0x4a,0x7b,0xb1,0xfe,0xe0,0x8f,0x57,0x12,0x42,0x42,0x50,0x51,0xc1,0xcc,0xd1,0x7f,0x9a,0xca,0xe0,0x1f,0x50,0x57};
    std::uint8_t product7[32]      = {0xfe,0x21,0xc4,0x1a,0x1e,0x0a,0x3d,0xb5,0x62,0xbd,0x89,0xad,0x80,0x8b,0xb4,0x29,0xdc,0x5a,0x0e,0xbe,0x98,0xf4,0xc1,0xbf,0xf0,0x9b,0x43,0x08,0x52,0x69,0xa5,0xc1};
    mont256_from_bytes(&ctx, multiplier, multiplier7);
    mont256_from_bytes(&ctx, multiplicand, multiplicand7);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product7, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier8[32]   = {0x45,0x1a,0xbd,0x81,0xf1,0xd6,0x9e,0xd6,0x17,0xf5,0xe8,0x37,0xd7,0x08,0x20,0xfe,0x11,0x9a,0x72,0xd1,0x74,0xc9,0xdf,0x6a,0xcc,0x01,0x1c,0xdd,0x94,0x74,0x03,0x1b};
    std::uint8_t multiplicand8[32] = {0x4f,0x42,0x6d,0xcb,0xb3,0x94,0xfb,0x36,0xbb,0x2d,0x42,0x0f,0x0f,0x88,0x08,0x0b,0x10,0xa3,0xd6,0xb2,0xaa,0x05,0xe1,0x1a,0xb2,0x71,0x59,0x45,0x79,0x5e,0x82,0x29};
    std::uint8_t product8[32]      = {0xb1,0xd4,0x63,0x58,0x42,0xf6,0x6e,0x9f,0x04,0x20,0x6a,0x1f,0xf0,0xd8,0x1f,0x87,0x82,0xbc,0x96,0x40,0x38,0xe2,0x54,0x42,0xd0,0xd6,0x68,0x44,0x6d,0x21,0x70,0xee};
    mont256_from_bytes(&ctx, multiplier, multiplier8);
    mont256_from_bytes(&ctx, multiplicand, multiplicand8);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product8, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier9[32]   = {0xb7,0x74,0xeb,0x52,0x48,0xdb,0x40,0xaf,0x72,0x15,0x83,0x70,0xd2,0x69,0xa9,0xa5,0xae,0x65,0x8f,0x33,0xfe,0x3b,0x89,0x0b,0x93,0xf4,0x48,0xb3,0xa5,0xaa,0x3c,0x81};
    std::uint8_t multiplicand9[32] = {0x5a,0xff,0xb2,0x29,0x76,0x31,0xa9,0x92,0xf0,0xce,0x58,0x35,0x05,0xc6,0xaf,0x07,0x58,0xd5,0x56,0x3d,0xab,0x2c,0xd3,0x1e,0xe3,0x15,0x12,0x88,0x62,0xc3,0x3a,0x4f};
    std::uint8_t product9[32]      = {0x8f,0xe2,0xf7,0x8d,0x50,0xd0,0x4a,0x98,0xc0,0x38,0x57,0xe9,0xce,0xad,0xe9,0x77,0x7b,0x7a,0x1a,0x23,0x45,0x8a,0xf6,0x88,0x90,0xa1,0x67,0x4e,0xaf,0xd0,0x9d,0x89};
    mont256_from_bytes(&ctx, multiplier, multiplier9);
    mont256_from_bytes(&ctx, multiplicand, multiplicand9);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product9, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier10[32]   = {0x49,0x95,0x23,0x99,0xc4,0xaa,0xea,0xc1,0x37,0xdc,0x76,0xfb,0x0f,0x17,0xa3,0x00,0x7e,0x62,0xaa,0x0a,0x1d,0xf9,0xfd,0x78,0x9c,0x65,0x39,0x38,0x2b,0x05,0x37,0xe6};
    std::uint8_t multiplicand10[32] = {0x7f,0x1b,0x10,0x3c,0xdf,0x15,0x82,0xb0,0xea,0xb4,0x77,0xd2,0x64,0x15,0x47,0x9c,0x65,0xdc,0x9f,0x50,0x3f,0x63,0xaf,0x83,0xbd,0x05,0x61,0xe6,0x21,0x1c,0x70,0xcf};
    std::uint8_t product10[32]      = {0x2b,0xa6,0x9d,0x52,0x53,0xae,0x98,0xd5,0x07,0xd1,0xce,0x00,0x50,0xb6,0x26,0x41,0x83,0xc3,0xef,0xd1,0xad,0x02,0x8a,0x93,0x4a,0x92,0xe3,0x82,0x86,0x32,0x1d,0xb7};
    mont256_from_bytes(&ctx, multiplier, multiplier10);
    mont256_from_bytes(&ctx, multiplicand, multiplicand10);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product10, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier11[32]   = {0x23,0x0d,0x97,0x7e,0xe2,0x25,0x71,0x59,0x47,0x20,0x77,0x1f,0x8c,0xa8,0x18,0x11,0x66,0xd2,0x28,0x76,0x72,0xfd,0xf2,0x02,0x2a,0x96,0xfb,0x1a,0x14,0xa0,0xf9,0xe7};
    std::uint8_t multiplicand11[32] = {0xfc,0x89,0x1b,0x4a,0x6a,0x50,0xdf,0x4d,0xb4,0xd6,0x6a,0x3a,0x47,0x46,0x9a,0x4d,0x8c,0xdb,0x30,0x5f,0xdd,0x2e,0x16,0x09,0x6e,0x36,0xaa,0xb0,0xd1,0xbc,0x52,0xd9};
    std::uint8_t product11[32]      = {0xf2,0xdf,0xfd,0x10,0x1f,0xd4,0xf7,0x2d,0xaa,0x10,0xf2,0xa1,0xae,0xc1,0x73,0x4e,0xef,0xec,0x30,0xd6,0xe3,0x19,0xda,0x21,0x76,0xd2,0x7b,0xdd,0x76,0xba,0xfb,0x7c};
    mont256_from_bytes(&ctx, multiplier, multiplier11);
    mont256_from_bytes(&ctx, multiplicand, multiplicand11);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product11, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier12[32]   = {0x15,0x3e,0x7c,0x2a,0x26,0xa2,0xc0,0xbd,0x3b,0x12,0x87,0xff,0xf5,0x2d,0xdf,0x5d,0x61,0x64,0x99,0xc9,0xe2,0x5a,0x76,0x05,0xae,0xc6,0xf0,0x24,0x5b,0xd8,0x6d,0x40};
    std::uint8_t multiplicand12[32] = {0xd4,0xc2,0x8c,0x2e,0x7c,0x26,0x84,0x7f,0x03,0x16,0x90,0x9e,0x3b,0xbb,0xe9,0xea,0xa8,0x94,0x8c,0x89,0x3b,0x61,0x86,0x76,0x26,0xbb,0x7d,0xbd,0x2d,0x1c,0x9a,0xf0};
    std::uint8_t product12[32]      = {0xcc,0xf9,0x39,0x2f,0xf2,0x7c,0x05,0x6d,0xc0,0x53,0x7a,0x4f,0x23,0x41,0x64,0x21,0xb0,0xe1,0x75,0xab,0x31,0xbc,0xf9,0x93,0x5e,0xa8,0x80,0xd6,0xb0,0x59,0x6c,0x3c};
    mont256_from_bytes(&ctx, multiplier, multiplier12);
    mont256_from_bytes(&ctx, multiplicand, multiplicand12);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product12, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier13[32]   = {0x88,0xda,0xf4,0x01,0x6b,0x40,0x13,0xef,0x25,0x4b,0x0c,0x4e,0x01,0x0c,0x47,0x59,0x48,0x2c,0x9c,0xbc,0x43,0x43,0x5c,0xc5,0x2e,0xae,0x05,0xcf,0x96,0xd0,0xcc,0x5f};
    std::uint8_t multiplicand13[32] = {0xdb,0xf4,0xa8,0xb2,0xb0,0xc4,0x31,0x2d,0x20,0x20,0x36,0x26,0xf3,0xfe,0x39,0xc0,0x51,0x90,0x88,0xf5,0x90,0xfb,0xbd,0x11,0x9c,0x1c,0xaa,0xf7,0x5e,0x87,0x66,0xed};
    std::uint8_t product13[32]      = {0x13,0x3c,0xe5,0x11,0x50,0x81,0xc9,0x58,0x36,0xfa,0x28,0x47,0xda,0xfd,0x4f,0x97,0x7a,0xe2,0xe8,0x8d,0xb9,0x36,0xd4,0xdc,0x51,0xf6,0x8e,0x7c,0xd1,0x1c,0x03,0x2b};
    mont256_from_bytes(&ctx, multiplier, multiplier13);
    mont256_from_bytes(&ctx, multiplicand, multiplicand13);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product13, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier14[32]   = {0x74,0xe6,0x9a,0x5d,0x0d,0xd2,0x7a,0x65,0xbd,0x62,0x88,0x81,0xad,0x1b,0x72,0xdb,0xa7,0xab,0xe1,0xc2,0x9e,0x1a,0x8e,0xf4,0xf3,0x41,0xe0,0x7a,0x83,0xf7,0x3f,0x16};
    std::uint8_t multiplicand14[32] = {0x8f,0x2c,0x6e,0xc8,0xcc,0x41,0x69,0xa3,0xae,0x3a,0x2b,0x7f,0xdf,0xe0,0x18,0x93,0xf3,0xae,0xd0,0xb6,0xc7,0xac,0x14,0x91,0xde,0xf8,0x83,0x34,0xe6,0x47,0xcb,0x8f};
    std::uint8_t product14[32]      = {0xce,0x93,0x33,0x94,0xb8,0xac,0x2f,0xed,0x2d,0xa9,0xd6,0xb9,0xce,0x11,0x2a,0x50,0x91,0x3d,0x52,0x8b,0xa7,0x5a,0x70,0xac,0x16,0x13,0x20,0x3f,0x36,0x7f,0x92,0x45};
    mont256_from_bytes(&ctx, multiplier, multiplier14);
    mont256_from_bytes(&ctx, multiplicand, multiplicand14);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product14, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier15[32]   = {0x66,0x83,0x68,0x86,0xa2,0x60,0xcd,0x0b,0x7b,0x45,0x14,0x5c,0x1a,0x81,0x68,0x2c,0x64,0xe5,0x0c,0xad,0x66,0x23,0x7a,0x04,0x65,0xe7,0xe4,0x23,0x64,0x72,0xf1,0xa3};
    std::uint8_t multiplicand15[32] = {0x1c,0x24,0x42,0xf9,0x29,0x8c,0xb3,0xa5,0x70,0xcc,0xec,0x31,0x35,0x71,0x81,0x0a,0xfc,0x13,0x2d,0x0d,0x11,0x3d,0xb1,0x7d,0x30,0xcb,0xc9,0x7d,0x0f,0xef,0x79,0x28};
    std::uint8_t product15[32]      = {0x18,0x53,0x16,0xf3,0x9b,0x93,0xaf,0xe6,0xb3,0xd4,0x6f,0xdf,0xbb,0xa5,0x63,0x28,0xdd,0x30,0x83,0xc0,0x92,0x80,0xa7,0xe3,0xef,0x37,0xd4,0x13,0x75,0xa1,0x2b,0x0b};
    mont256_from_bytes(&ctx, multiplier, multiplier15);
    mont256_from_bytes(&ctx, multiplicand, multiplicand15);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product15, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier16[32]   = {0x89,0x5f,0xd7,0xb3,0x26,0xb9,0x4c,0x7f,0x91,0x18,0xbb,0x16,0x00,0x0f,0x49,0xc8,0x1a,0x35,0x8c,0xa0,0x0d,0x75,0x98,0x5d,0x99,0xc9,0x43,0x09,0x57,0x0d,0xc1,0x95};
    std::uint8_t multiplicand16[32] = {0x35,0x3c,0x63,0x1c,0xdf,0xd4,0x3f,0x37,0x12,0x00,0x33,0x9d,0x06,0x87,0x39,0xfa,0x9d,0x1d,0xe2,0xa0,0x5d,0x15,0x8a,0x2f,0xf2,0xee,0x4e,0x45,0x19,0xf9,0x91,0x9c};
    std::uint8_t product16[32]      = {0x16,0xd3,0xa7,0x19,0x5d,0xf1,0x54,0x0b,0x8f,0x2c,0x70,0x61,0x6a,0x9e,0xb5,0xdf,0xac,0xfc,0x82,0x86,0x39,0x94,0x68,0x62,0x5d,0x33,0x55,0x07,0xae,0xde,0xdc,0xcc};
    mont256_from_bytes(&ctx, multiplier, multiplier16);
    mont256_from_bytes(&ctx, multiplicand, multiplicand16);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product16, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier17[32]   = {0x9a,0x2e,0xf8,0x0f,0x58,0xee,0x85,0x71,0xf4,0x99,0x8d,0x7c,0x40,0x93,0xf6,0xde,0xa2,0x68,0xaa,0x87,0x26,0x07,0x67,0x9d,0x60,0x50,0x91,0x4a,0x9d,0x33,0xa0,0x1c};
    std::uint8_t multiplicand17[32] = {0xfa,0x52,0x9b,0xa3,0xfe,0x3b,0xfa,0xda,0x7c,0xf2,0x07,0x24,0xd9,0x53,0xee,0x26,0x1d,0x87,0xce,0xc3,0x1f,0x72,0x96,0xab,0x79,0x61,0xfd,0x92,0x5d,0x39,0xd0,0xa8};
    std::uint8_t product17[32]      = {0xb1,0x3e,0xbe,0xee,0x7b,0xf7,0x9b,0xb3,0xdf,0x93,0xca,0x05,0x6a,0x79,0x39,0xa5,0xd0,0xb3,0x74,0xaf,0x1c,0x81,0x25,0x79,0x26,0xbc,0xbe,0x58,0x7f,0xae,0xb3,0x95};
    mont256_from_bytes(&ctx, multiplier, multiplier17);
    mont256_from_bytes(&ctx, multiplicand, multiplicand17);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product17, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier18[32]   = {0xbf,0xea,0xa1,0x55,0x1a,0x28,0xf7,0xb3,0x24,0xe4,0xe2,0x5a,0x15,0xfc,0x89,0x9e,0x4f,0xd5,0x8d,0xbe,0x7b,0xdc,0x96,0x8b,0x7a,0xfb,0x2c,0x68,0x77,0x4b,0x15,0xd7};
    std::uint8_t multiplicand18[32] = {0x84,0x2e,0x7f,0xc2,0x29,0x54,0x0a,0x6e,0xb1,0x2a,0xa1,0xf6,0xd4,0x2f,0xdd,0xbb,0x7a,0x86,0xf7,0xa2,0x43,0xc7,0x1b,0x9a,0xbd,0x87,0xa8,0x65,0x57,0xb6,0xfb,0x7e};
    std::uint8_t product18[32]      = {0x18,0xfa,0xb7,0x21,0x92,0xbb,0xfe,0xe0,0xe9,0x12,0x1b,0xb2,0xe5,0xc1,0x76,0x8e,0x5c,0x47,0xd3,0x2a,0x2e,0xb2,0xe9,0x7d,0xa4,0x17,0xb6,0x8c,0x6f,0x9d,0xfd,0x5d};
    mont256_from_bytes(&ctx, multiplier, multiplier18);
    mont256_from_bytes(&ctx, multiplicand, multiplicand18);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product18, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier19[32]   = {0xb0,0xa8,0x44,0xe5,0x25,0x87,0xbe,0x6b,0x5c,0x9b,0xcf,0x35,0x87,0x3b,0xe0,0x78,0xf3,0xb7,0xa5,0x0d,0xf3,0x73,0xca,0x53,0x34,0x88,0xf8,0x76,0x05,0xe9,0x99,0xf3};
    std::uint8_t multiplicand19[32] = {0xa4,0x96,0x36,0xa2,0xfa,0x7f,0x0e,0xab,0x4c,0x4f,0x9b,0x06,0x87,0x32,0x2e,0x25,0xc2,0x15,0xa8,0x2a,0x06,0xec,0x41,0xad,0xea,0x05,0x75,0x43,0x8b,0x0d,0x59,0x0b};
    std::uint8_t product19[32]      = {0x1f,0xc0,0x3e,0x32,0x5b,0xd7,0x8a,0xe9,0x88,0x09,0xfb,0x1b,0xcd,0x17,0x29,0xd7,0x76,0x75,0x57,0x25,0xb9,0x59,0x25,0x4c,0x75,0x1a,0xe8,0xe4,0xcb,0xa0,0x8e,0xf7};
    mont256_from_bytes(&ctx, multiplier, multiplier19);
    mont256_from_bytes(&ctx, multiplicand, multiplicand19);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product19, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier20[32]   = {0xe8,0x83,0xa1,0xd4,0x5d,0xe0,0x09,0x97,0x84,0xb5,0xa8,0x18,0x42,0xd8,0x72,0x08,0xd8,0x6f,0x40,0xf6,0xb2,0x39,0xf3,0xc7,0x17,0x4c,0x77,0xa2,0xdd,0x02,0xde,0x92};
    std::uint8_t multiplicand20[32] = {0x80,0xb0,0xc0,0x8b,0xc7,0x70,0x24,0x20,0x8a,0xa4,0x24,0x8c,0x88,0x57,0xf9,0xa4,0x39,0x08,0xf2,0x27,0xc5,0x9d,0xb9,0x16,0x5b,0x0e,0xe7,0x6f,0x2a,0xc3,0x44,0x46};
    std::uint8_t product20[32]      = {0xaa,0x4e,0xe4,0x09,0xd7,0x51,0xa1,0x6e,0x48,0xac,0x7f,0xd7,0x6d,0x2c,0xe7,0xfe,0x97,0xb3,0xac,0xef,0x6a,0xd8,0xec,0x6b,0xca,0x61,0x72,0x7d,0xd3,0x4c,0xb2,0xf7};
    mont256_from_bytes(&ctx, multiplier, multiplier20);
    mont256_from_bytes(&ctx, multiplicand, multiplicand20);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product20, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier21[32]   = {0xc2,0x21,0x6b,0x02,0xfc,0x24,0x1d,0x0b,0xc9,0xd4,0x88,0xb1,0xcf,0xbf,0x33,0x60,0x9c,0xfc,0x86,0x52,0x39,0x19,0x42,0x42,0xa2,0xed,0xdb,0xbd,0x54,0x64,0xec,0xc2};
    std::uint8_t multiplicand21[32] = {0xcd,0xa6,0xc6,0xfd,0xbd,0x68,0x51,0x67,0x66,0x93,0x40,0x36,0xd1,0x7e,0x44,0x97,0x3d,0x48,0x82,0xa5,0xce,0x5b,0x2a,0x92,0x31,0xf5,0x17,0x07,0xda,0x45,0xe1,0x8a};
    std::uint8_t product21[32]      = {0x35,0xed,0xf6,0xbe,0xd6,0xa7,0x2d,0xbf,0x79,0xe3,0x1f,0x69,0xe5,0x5a,0xde,0x48,0xb9,0x5d,0x00,0xeb,0xe6,0x13,0x90,0x6c,0x49,0x7e,0xa5,0xa0,0x39,0x10,0xa1,0xd2};
    mont256_from_bytes(&ctx, multiplier, multiplier21);
    mont256_from_bytes(&ctx, multiplicand, multiplicand21);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product21, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier22[32]   = {0xfd,0x56,0xa9,0x26,0x07,0x6b,0x3e,0x36,0xbb,0x23,0x13,0xf5,0x5b,0x06,0x25,0x8e,0x7e,0x26,0xf3,0x6a,0x84,0x83,0xf8,0xb8,0x33,0x2d,0xd3,0x31,0x3a,0x0b,0x99,0x65};
    std::uint8_t multiplicand22[32] = {0x9a,0xea,0x64,0x29,0xb1,0x49,0x1e,0x24,0x31,0x92,0xb7,0x04,0x42,0x59,0x40,0x52,0x78,0xe4,0xb9,0x8d,0x47,0x87,0xf9,0x3b,0xca,0x44,0xeb,0x86,0x07,0x26,0xe2,0x5c};
    std::uint8_t product22[32]      = {0x9a,0x54,0x51,0xc9,0x4a,0xa6,0x60,0x99,0xbe,0xa9,0x52,0x8a,0xef,0xfa,0xe8,0x0b,0x58,0x54,0xcf,0xb4,0xec,0xb4,0x38,0xf4,0x4c,0x52,0x4e,0x3b,0x45,0x9b,0x79,0x17};
    mont256_from_bytes(&ctx, multiplier, multiplier22);
    mont256_from_bytes(&ctx, multiplicand, multiplicand22);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product22, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier23[32]   = {0x59,0x7a,0x1e,0xcf,0xfc,0xf0,0x0f,0xec,0xb9,0x1e,0xe9,0xe5,0xef,0xe0,0x9f,0x07,0xce,0xfe,0x2a,0x1f,0x72,0x7d,0x83,0x49,0x58,0x22,0xcb,0x77,0xf4,0xde,0x2c,0x08};
    std::uint8_t multiplicand23[32] = {0x78,0x57,0x29,0x76,0x3a,0x12,0x91,0x7c,0x1a,0x26,0xf8,0x89,0x38,0x70,0x38,0x00,0x14,0x9e,0x25,0x9b,0x5d,0x58,0xc7,0x05,0xf9,0x79,0xd0,0x4a,0xf4,0x7a,0xeb,0xdd};
    std::uint8_t product23[32]      = {0x9f,0x36,0xc7,0xd2,0xaa,0x85,0x71,0x50,0x2f,0x32,0x47,0xf8,0x04,0xc5,0x69,0x75,0x45,0xe9,0xf2,0xac,0x43,0x6a,0x02,0x58,0xa5,0x4f,0xef,0x19,0xa3,0x56,0x8f,0xbf};
    mont256_from_bytes(&ctx, multiplier, multiplier23);
    mont256_from_bytes(&ctx, multiplicand, multiplicand23);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product23, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier24[32]   = {0x9c,0x3a,0x23,0xcd,0xe6,0x7a,0x9b,0x75,0xfc,0x39,0x47,0x24,0x9f,0xc2,0xd0,0xa1,0x7b,0x8f,0x2a,0xb5,0x34,0x51,0xd0,0x13,0x56,0x75,0xf6,0xad,0x32,0x5b,0x55,0xdd};
    std::uint8_t multiplicand24[32] = {0xa4,0xa4,0x5e,0xff,0xcc,0xb5,0x73,0xd9,0x58,0x10,0xd6,0x0e,0xa7,0x29,0x91,0xb9,0xe8,0xc1,0x47,0x43,0x7a,0xbe,0xc5,0x39,0x00,0x7d,0x10,0x34,0xd7,0x26,0xc8,0x6b};
    std::uint8_t product24[32]      = {0xac,0x6a,0xf5,0x11,0xc6,0xea,0xbf,0xd3,0x6a,0x9a,0x2a,0x28,0x7c,0x13,0x94,0x1f,0x40,0xf2,0x67,0x6d,0xf8,0x1b,0x38,0xfe,0x5c,0x68,0x9a,0x4a,0x27,0xfd,0x42,0x32};
    mont256_from_bytes(&ctx, multiplier, multiplier24);
    mont256_from_bytes(&ctx, multiplicand, multiplicand24);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product24, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier25[32]   = {0xb6,0x24,0x67,0x71,0xc8,0x45,0x00,0x70,0x63,0x77,0x14,0x07,0xe8,0xe7,0x27,0x89,0x1e,0xb2,0x01,0x09,0xa9,0x1c,0x24,0x39,0xd5,0xab,0x8b,0x4d,0x15,0xb4,0x0a,0xeb};
    std::uint8_t multiplicand25[32] = {0xa2,0xc6,0x8e,0x45,0xca,0x04,0xc7,0x9f,0x6f,0x15,0xb6,0xad,0x2d,0xb3,0x99,0x7f,0xe3,0x96,0x39,0xbe,0x7a,0x60,0x5a,0x91,0x33,0x06,0x98,0xa1,0xc0,0x09,0x34,0x92};
    std::uint8_t product25[32]      = {0x66,0x1c,0x82,0x0a,0x05,0xf4,0x46,0x49,0x1e,0x72,0x26,0xa9,0x19,0xb1,0x80,0xc5,0xcd,0xc6,0xb9,0x21,0xb5,0x62,0xd4,0xb2,0x0a,0x8d,0x3d,0xdc,0x1b,0xf5,0x12,0xe5};
    mont256_from_bytes(&ctx, multiplier, multiplier25);
    mont256_from_bytes(&ctx, multiplicand, multiplicand25);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product25, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier26[32]   = {0x76,0x91,0xb0,0x6f,0x65,0x55,0xab,0xfe,0xb8,0xc9,0x81,0x7a,0xf8,0xbe,0x88,0x31,0xf2,0x37,0xe4,0x5a,0xcd,0x02,0xc5,0xe1,0x16,0x35,0x3d,0x03,0x55,0x1f,0xd8,0xf9};
    std::uint8_t multiplicand26[32] = {0xfe,0x3c,0x9c,0x8f,0x2b,0x85,0x5c,0x1f,0x28,0xaa,0xca,0x51,0xb9,0x8c,0x67,0xc2,0x15,0xbd,0x44,0x8f,0xf2,0x61,0x49,0xed,0xbe,0x4c,0x5c,0xe6,0x66,0xc1,0x49,0x4e};
    std::uint8_t product26[32]      = {0xb2,0x9b,0xc4,0x5d,0xcd,0xe3,0xec,0x96,0xb4,0xca,0xb8,0xae,0x6d,0x16,0x68,0xad,0x7c,0x69,0x48,0x5d,0xfd,0x56,0xd8,0xbb,0x8e,0x64,0xdf,0x56,0xea,0x83,0x0b,0xf7};
    mont256_from_bytes(&ctx, multiplier, multiplier26);
    mont256_from_bytes(&ctx, multiplicand, multiplicand26);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product26, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier27[32]   = {0xa7,0xe6,0x52,0x9b,0xce,0x76,0xe9,0xf4,0x77,0x21,0x6e,0x9e,0xe7,0xa4,0x63,0x09,0x97,0x3f,0x79,0x86,0x26,0xb1,0xcf,0xfc,0x07,0x0d,0x71,0x09,0x20,0x85,0x96,0x34};
    std::uint8_t multiplicand27[32] = {0xef,0xfd,0xde,0xea,0xa8,0x42,0xbc,0x19,0x79,0x6f,0x74,0xad,0xfa,0xf5,0x54,0x96,0x98,0x8a,0xf3,0xfb,0xd3,0x96,0x30,0xd6,0x9c,0x90,0x11,0xef,0x25,0x6b,0xad,0xf9};
    std::uint8_t product27[32]      = {0x4b,0x00,0x18,0x82,0xcd,0xcd,0x79,0x35,0x7b,0x3f,0x60,0x5c,0xbe,0x55,0xde,0x01,0x09,0xba,0xb2,0x33,0x83,0x3a,0x5e,0xe8,0xe8,0x54,0xd2,0x9e,0x02,0x28,0x29,0xc9};
    mont256_from_bytes(&ctx, multiplier, multiplier27);
    mont256_from_bytes(&ctx, multiplicand, multiplicand27);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product27, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier28[32]   = {0xcc,0xa2,0xa9,0x2b,0x03,0xa5,0x6c,0xc1,0x05,0x7a,0x40,0xb2,0x21,0x88,0x28,0x7e,0x8c,0x5c,0x71,0x5f,0x8c,0x74,0xfc,0x1e,0x27,0xe9,0xe0,0x6f,0x59,0xb4,0x4e,0x92};
    std::uint8_t multiplicand28[32] = {0x23,0xa5,0xef,0x88,0xef,0x02,0x09,0x0b,0xbf,0xde,0xfc,0x15,0x86,0xce,0x03,0xf9,0x1a,0x4f,0x44,0xf9,0xa6,0x51,0x14,0x45,0xb9,0xf3,0x63,0x5c,0xf8,0x8c,0x42,0x2b};
    std::uint8_t product28[32]      = {0x2b,0xf4,0x78,0x7e,0x88,0xe4,0xb4,0x6f,0x8c,0x75,0x53,0x3b,0x29,0xf3,0xcd,0x6a,0xa1,0x12,0xaf,0x5e,0x76,0x85,0x15,0xe9,0x6d,0xb7,0xa2,0x89,0x3f,0x76,0x84,0x8b};
    mont256_from_bytes(&ctx, multiplier, multiplier28);
    mont256_from_bytes(&ctx, multiplicand, multiplicand28);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product28, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier29[32]   = {0x07,0x2a,0x98,0xd2,0x36,0x06,0xde,0xfc,0xdf,0xb8,0x5c,0x0d,0xd3,0x7e,0xe9,0x15,0x31,0xde,0xc4,0xf4,0xdf,0x2a,0x8b,0x79,0xfc,0x8e,0x80,0xb3,0x6f,0x0e,0x22,0x89};
    std::uint8_t multiplicand29[32] = {0x53,0x74,0x09,0x02,0x96,0x20,0xbf,0x0d,0xc3,0x80,0x84,0xa0,0x3d,0x93,0xfd,0x4c,0x80,0x4c,0x25,0xd6,0x4a,0xff,0xdc,0xd1,0x36,0x78,0xbc,0x8d,0x40,0x78,0x3f,0x0a};
    std::uint8_t product29[32]      = {0x3c,0x7e,0xba,0x56,0x3b,0xef,0x9e,0xee,0x50,0xa0,0xfb,0xf3,0xa5,0xdb,0xbb,0x9c,0x0e,0xab,0xfc,0xef,0x9d,0x5d,0x89,0xa1,0xfb,0xca,0xd2,0x92,0x8e,0x5f,0x4c,0xc6};
    mont256_from_bytes(&ctx, multiplier, multiplier29);
    mont256_from_bytes(&ctx, multiplicand, multiplicand29);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product29, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier30[32]   = {0xbd,0x6b,0x88,0x1a,0xe8,0xf6,0xe0,0xbd,0x0f,0x97,0x70,0x44,0x21,0x8e,0x0b,0x7b,0xd5,0x8d,0xcd,0xb4,0x6b,0x44,0x68,0x06,0x8b,0x5a,0xb3,0xee,0x42,0x65,0xbb,0x31};
    std::uint8_t multiplicand30[32] = {0x84,0x4a,0x70,0x34,0xe7,0x7f,0xfe,0x48,0xd0,0xa6,0xec,0x17,0x95,0x56,0x58,0x5e,0xa9,0x97,0xf3,0x51,0x75,0x4a,0x09,0xcd,0xe5,0xcf,0xed,0xfa,0x5a,0x91,0x96,0xf0};
    std::uint8_t product30[32]      = {0x15,0x29,0x23,0x8f,0x17,0x28,0x0c,0xdc,0xa2,0x5d,0x62,0x14,0xc0,0x82,0x9a,0x78,0xa2,0x99,0xa6,0x5c,0x53,0xfc,0xcb,0xd5,0x39,0xce,0x2e,0x0e,0x93,0x93,0x88,0x37};
    mont256_from_bytes(&ctx, multiplier, multiplier30);
    mont256_from_bytes(&ctx, multiplicand, multiplicand30);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product30, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier31[32]   = {0x26,0xde,0xbf,0xdb,0x88,0x25,0xae,0x56,0x21,0x79,0xb3,0x7d,0x80,0x6c,0x10,0xb5,0xe0,0xcf,0xab,0x4c,0xea,0xef,0xc4,0xd2,0xd3,0xbf,0x6d,0x01,0x6b,0xae,0x4b,0x5b};
    std::uint8_t multiplicand31[32] = {0x9b,0xca,0x3c,0xb7,0x2e,0xe0,0x28,0x9d,0xc6,0xc9,0x1b,0x92,0x70,0xac,0x06,0xac,0xdf,0x70,0x30,0x17,0x04,0xc9,0xd7,0x8d,0x82,0xb3,0x35,0x99,0x86,0x04,0x87,0x19};
    std::uint8_t product31[32]      = {0xe4,0x2f,0xd8,0xe4,0xa2,0x2f,0xb9,0x7e,0xe0,0xa5,0xf8,0x0c,0x75,0xc0,0xd6,0xc7,0xe2,0xb8,0xd3,0x71,0x3d,0x57,0x0f,0x58,0x2d,0xe5,0xc0,0x5c,0x06,0x2a,0xf0,0x4c};
    mont256_from_bytes(&ctx, multiplier, multiplier31);
    mont256_from_bytes(&ctx, multiplicand, multiplicand31);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product31, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier32[32]   = {0x9e,0x7d,0x6b,0x37,0x79,0x36,0xd5,0x36,0x24,0x3d,0x35,0x70,0x2c,0x1e,0xea,0x1f,0x26,0x59,0x74,0xa7,0xcc,0x96,0x6f,0x46,0xc6,0xaa,0x7d,0x55,0x01,0x01,0xb8,0x11};
    std::uint8_t multiplicand32[32] = {0x87,0xdd,0xae,0xb7,0x84,0xb2,0x80,0x54,0xae,0xad,0x44,0xb0,0x53,0x73,0x90,0xe5,0x0f,0xcf,0x31,0xca,0x8e,0x75,0x2f,0xdf,0x1e,0xce,0x61,0x5d,0xb9,0xa6,0x44,0x2e};
    std::uint8_t product32[32]      = {0x80,0xae,0x3e,0xc5,0xe1,0xad,0x96,0x48,0x3d,0xd2,0x50,0xc2,0x77,0x4e,0x8a,0xa8,0x95,0xaf,0x92,0x0a,0x6a,0x3b,0xb9,0x8c,0x9a,0x2d,0x72,0x6a,0xb6,0x63,0xda,0x7b};
    mont256_from_bytes(&ctx, multiplier, multiplier32);
    mont256_from_bytes(&ctx, multiplicand, multiplicand32);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product32, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier33[32]   = {0x0e,0x8b,0xec,0x94,0x8f,0x6f,0x91,0x5f,0xe2,0x1b,0x37,0xca,0x1b,0x29,0xfc,0x99,0xc6,0xc8,0x0e,0x2b,0xc8,0xc6,0x14,0xb2,0x7b,0x84,0x44,0xd1,0x8e,0x31,0x70,0x41};
    std::uint8_t multiplicand33[32] = {0x73,0xc1,0xcd,0x2c,0x81,0xf9,0x8b,0x52,0x19,0x05,0xd5,0x91,0xc5,0xb2,0xe7,0x5a,0x0a,0xcd,0x8b,0xe1,0x46,0xe4,0x09,0x90,0x30,0xf9,0x70,0x58,0x3f,0x9d,0x52,0xf9};
    std::uint8_t product33[32]      = {0x81,0x61,0xdc,0x50,0x74,0xca,0x99,0xa3,0xdc,0x5f,0x09,0xac,0x1d,0xed,0x9c,0x2a,0x39,0x1c,0x9a,0x67,0xd1,0x91,0xb1,0x02,0xdd,0xff,0x5c,0x51,0xa2,0x7d,0x9e,0x22};
    mont256_from_bytes(&ctx, multiplier, multiplier33);
    mont256_from_bytes(&ctx, multiplicand, multiplicand33);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product33, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier34[32]   = {0x53,0x5b,0x6a,0x43,0x71,0x78,0xba,0x0a,0x10,0x38,0xf0,0xb5,0xe9,0x98,0xd0,0xee,0xe4,0xdd,0xf9,0xb9,0xc2,0x8e,0xe9,0x07,0x07,0x22,0x35,0xc2,0x8f,0xcd,0x7f,0x40};
    std::uint8_t multiplicand34[32] = {0x46,0xf5,0xa1,0xb4,0xb1,0x56,0xd1,0xad,0x33,0x0c,0x16,0xa3,0x83,0x1d,0x03,0xbf,0x9b,0x2b,0xd6,0xc0,0x81,0x6b,0xee,0x06,0xf9,0x2e,0x23,0x39,0x9c,0xce,0xa0,0x98};
    std::uint8_t product34[32]      = {0x17,0x00,0x26,0x7f,0x4e,0x57,0x23,0x73,0x04,0xc2,0xa3,0xf1,0x6e,0x95,0xb9,0x05,0x7b,0xa9,0x04,0x68,0x19,0x93,0x5b,0xb6,0x9e,0xf0,0xe4,0x69,0x07,0x5d,0xf9,0xa9};
    mont256_from_bytes(&ctx, multiplier, multiplier34);
    mont256_from_bytes(&ctx, multiplicand, multiplicand34);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product34, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier35[32]   = {0x3f,0x66,0x5e,0xde,0xf1,0x06,0x37,0xce,0x81,0xfc,0x06,0x9e,0x7a,0x60,0x96,0x83,0xce,0xaf,0x49,0x15,0x88,0x85,0x64,0xe8,0x82,0x16,0x85,0x8f,0x73,0xcc,0xef,0x03};
    std::uint8_t multiplicand35[32] = {0xec,0x3b,0x96,0x05,0x42,0x74,0xa3,0xeb,0xed,0x84,0xe9,0x1e,0xf1,0x32,0xbf,0x2d,0xe0,0x40,0x01,0x5c,0xe0,0x64,0xa1,0x14,0x85,0xf1,0x11,0x5b,0xb2,0xff,0xf1,0x7b};
    std::uint8_t product35[32]      = {0x3b,0x43,0x1b,0xb1,0xfb,0xdc,0x0e,0x39,0x08,0xca,0xb0,0xde,0xdd,0x03,0xbd,0x04,0xe5,0x14,0x4d,0xd8,0x57,0x8d,0x38,0x7c,0x23,0x3f,0x90,0xe6,0xb9,0x75,0xbb,0xd4};
    mont256_from_bytes(&ctx, multiplier, multiplier35);
    mont256_from_bytes(&ctx, multiplicand, multiplicand35);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product35, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier36[32]   = {0x6a,0xa8,0xb9,0xe0,0x23,0x1b,0x3e,0x14,0x72,0x91,0x35,0xbd,0xd7,0x0a,0x39,0xd1,0x33,0xdc,0xd7,0x7f,0xf1,0x79,0xf2,0xd2,0xe4,0x8b,0x96,0x62,0x8f,0x3c,0x4b,0xe3};
    std::uint8_t multiplicand36[32] = {0x6d,0xa7,0x9a,0x87,0x3d,0x9a,0x80,0x79,0xab,0xd0,0xd7,0xfb,0x12,0x92,0x61,0x85,0x50,0xe4,0x0d,0x54,0x71,0x2e,0xa6,0xb3,0x64,0x71,0xfd,0xe4,0x1f,0x22,0x9d,0xd0};
    std::uint8_t product36[32]      = {0xb2,0x51,0xb0,0xe5,0x28,0xa7,0x1d,0xb8,0x7d,0x65,0x62,0xc5,0xfc,0xbb,0x54,0x45,0x48,0x6d,0x70,0x99,0x86,0x52,0xe8,0x64,0xa1,0x7c,0xac,0x98,0x7e,0xf5,0x62,0x84};
    mont256_from_bytes(&ctx, multiplier, multiplier36);
    mont256_from_bytes(&ctx, multiplicand, multiplicand36);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product36, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier37[32]   = {0xc6,0xe5,0x0d,0xf2,0xe5,0xa3,0x86,0x3e,0x1f,0x52,0x52,0x65,0xc8,0xb0,0x07,0xee,0x4d,0x82,0xfe,0xac,0xab,0x62,0x86,0xcd,0x36,0x72,0xd6,0xae,0x12,0xb8,0x0a,0xed};
    std::uint8_t multiplicand37[32] = {0x40,0xcb,0xac,0xd0,0x24,0x9a,0x45,0x84,0x5d,0xbe,0x30,0x23,0xa9,0x06,0x92,0x2f,0xa4,0xb9,0xa9,0xc4,0xb7,0x53,0xa1,0xee,0xf0,0x83,0x60,0x85,0x27,0x89,0xd0,0x59};
    std::uint8_t product37[32]      = {0xfb,0x04,0x4d,0xd4,0xff,0x99,0x0b,0xdc,0x3b,0xa2,0x4f,0x0a,0xc2,0x2f,0x33,0x63,0x70,0x3c,0xbd,0x1a,0xd5,0x21,0xce,0x99,0xaf,0x8d,0xd9,0x20,0x49,0x62,0x67,0x5b};
    mont256_from_bytes(&ctx, multiplier, multiplier37);
    mont256_from_bytes(&ctx, multiplicand, multiplicand37);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product37, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier38[32]   = {0x18,0x18,0x9a,0xf4,0xf3,0xd7,0x4f,0x82,0xbf,0x26,0x8e,0xa0,0x38,0x36,0xe8,0x65,0x77,0xbd,0x89,0x1f,0xf7,0xb1,0x03,0xdf,0x23,0x23,0x1e,0x1e,0xe2,0x01,0x55,0x22};
    std::uint8_t multiplicand38[32] = {0x39,0x45,0x33,0x6b,0xd5,0x1b,0x18,0x15,0xaa,0xf7,0x19,0xf3,0xfd,0x68,0x37,0x3b,0x29,0xac,0xf1,0xa5,0x7c,0xbd,0x1f,0x5a,0xe2,0x8a,0xf6,0x04,0x65,0xf4,0x29,0x86};
    std::uint8_t product38[32]      = {0x4f,0xd8,0xa6,0x75,0xf0,0xb2,0xfc,0xcc,0x25,0xf3,0xd0,0x52,0xe8,0x56,0x17,0x11,0x1b,0x2c,0x91,0x34,0xde,0xea,0xdc,0x9c,0xd3,0x9f,0x7c,0xdf,0x92,0x8a,0x8a,0xb1};
    mont256_from_bytes(&ctx, multiplier, multiplier38);
    mont256_from_bytes(&ctx, multiplicand, multiplicand38);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product38, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier39[32]   = {0x6b,0xd8,0xc6,0x76,0x56,0xd0,0x50,0xcd,0x67,0x60,0x13,0x67,0x83,0xfe,0xb1,0x7b,0xfe,0x7b,0x8a,0xe4,0x6e,0x78,0x36,0xa4,0xb4,0xd1,0x9e,0xc1,0x29,0x55,0xd6,0xf0};
    std::uint8_t multiplicand39[32] = {0x56,0x85,0xd6,0x24,0x04,0xfc,0xd5,0x55,0x5d,0xaf,0x10,0x6d,0xb8,0xde,0xe0,0x81,0x17,0x9a,0x07,0x1e,0x51,0x8a,0xe4,0x52,0x5b,0x4b,0x1b,0x75,0x32,0x1c,0x52,0x96};
    std::uint8_t product39[32]      = {0xf6,0xc7,0x56,0xee,0x5a,0xae,0x62,0x57,0xbd,0xf7,0x4a,0xe5,0x43,0xd0,0x6c,0x6a,0xa5,0x89,0xb0,0xb0,0xd1,0x3e,0x82,0x9d,0x77,0x8f,0xc1,0xdb,0x90,0xa9,0x50,0x87};
    mont256_from_bytes(&ctx, multiplier, multiplier39);
    mont256_from_bytes(&ctx, multiplicand, multiplicand39);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product39, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier40[32]   = {0x84,0x76,0x8b,0x8c,0x54,0xdd,0x0b,0xa5,0x62,0x64,0x67,0xba,0x04,0xa1,0x05,0x47,0xb4,0x01,0xba,0x85,0x70,0xc1,0xdc,0xa1,0x75,0x6b,0x72,0x89,0x8d,0xd6,0x3c,0xb9};
    std::uint8_t multiplicand40[32] = {0xeb,0x25,0xf8,0xa1,0xfc,0x2e,0x6a,0x59,0x1c,0xe3,0xbc,0x0c,0x10,0x75,0x5c,0x97,0xf5,0xf5,0x54,0xed,0x83,0x23,0x9e,0xf5,0x4b,0xa2,0xe1,0x61,0x9f,0xb9,0xaf,0x50};
    std::uint8_t product40[32]      = {0x05,0x4f,0x6f,0x54,0xe0,0xed,0x72,0x7c,0xce,0x1f,0x29,0x49,0x46,0xf8,0x12,0x46,0xa5,0xff,0x95,0xc8,0x1d,0x40,0x0e,0x4f,0xa5,0x44,0x85,0xf2,0xee,0x3e,0x42,0xc9};
    mont256_from_bytes(&ctx, multiplier, multiplier40);
    mont256_from_bytes(&ctx, multiplicand, multiplicand40);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product40, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier41[32]   = {0x45,0x9c,0x94,0x5c,0x43,0xfc,0x05,0x27,0x15,0x85,0x0a,0x03,0x1a,0xd2,0xd5,0xf1,0xe0,0x5b,0x3e,0x13,0xf8,0xc1,0x10,0xfb,0x3a,0x82,0x81,0x59,0xc9,0xd2,0x29,0x50};
    std::uint8_t multiplicand41[32] = {0xd1,0xdc,0xec,0x53,0x21,0x2a,0x8d,0x9b,0xc1,0x7a,0x92,0x62,0x45,0x3b,0xf4,0x91,0x2e,0x7a,0x26,0xe9,0xc7,0x6c,0x60,0x3f,0xe7,0xe8,0xf9,0xf6,0x0a,0x22,0x73,0x85};
    std::uint8_t product41[32]      = {0xc2,0x9f,0xfe,0x95,0x7d,0x90,0xff,0x64,0x25,0xe4,0xc4,0x17,0xd5,0x41,0xf3,0xf2,0xbd,0x24,0xe7,0x95,0x76,0xfb,0x5f,0x22,0x88,0x6d,0xc4,0x75,0xb2,0x0d,0x20,0xec};
    mont256_from_bytes(&ctx, multiplier, multiplier41);
    mont256_from_bytes(&ctx, multiplicand, multiplicand41);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product41, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier42[32]   = {0x67,0xec,0x32,0x6a,0x42,0x34,0x33,0x54,0xf2,0x2d,0x28,0x82,0xd1,0xa8,0x9b,0x37,0xad,0x0c,0x9b,0xb6,0xe9,0x52,0x6a,0x69,0xd9,0x7e,0x96,0x7b,0x6c,0x18,0xd9,0x82};
    std::uint8_t multiplicand42[32] = {0x53,0xb9,0x73,0x77,0xb3,0x4e,0x8e,0xce,0x7e,0x9e,0xe5,0x1d,0x92,0x12,0x82,0x4c,0x83,0xc8,0xcb,0x28,0xeb,0x4e,0xd2,0xe3,0x89,0x5e,0x8b,0x6b,0x26,0x3c,0xfa,0x5e};
    std::uint8_t product42[32]      = {0x6c,0x13,0x07,0xdb,0x92,0x4e,0x8d,0xda,0x62,0x25,0x6d,0xc3,0xaf,0xc4,0x30,0xb9,0x02,0x31,0x08,0xb4,0x4d,0x1a,0x80,0x42,0xa8,0x71,0xcf,0x2f,0xbf,0xc0,0x91,0x4d};
    mont256_from_bytes(&ctx, multiplier, multiplier42);
    mont256_from_bytes(&ctx, multiplicand, multiplicand42);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product42, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier43[32]   = {0xe5,0x31,0x69,0x60,0x6c,0xe1,0x93,0xc2,0x2e,0xef,0xa2,0x79,0xb0,0x2e,0x3d,0x8d,0xcc,0xb1,0xc5,0x1d,0x0e,0xba,0x0e,0xa8,0x47,0x70,0xa0,0x87,0x16,0xe6,0xfe,0xc3};
    std::uint8_t multiplicand43[32] = {0x42,0xb3,0x87,0x55,0xcd,0x37,0x88,0x0e,0x16,0xac,0x41,0x91,0xa2,0x6a,0xa0,0xae,0x04,0x4f,0x15,0x74,0xf0,0x37,0xaf,0xc6,0x44,0xd8,0x2a,0x53,0x12,0x89,0xba,0xfa};
    std::uint8_t product43[32]      = {0xff,0x57,0x51,0xca,0x5f,0xaf,0x1e,0x1a,0xb2,0x80,0x00,0xc6,0xa9,0x6e,0xa8,0x7f,0x1a,0xc0,0xe1,0xdf,0x3f,0x8b,0x21,0xdb,0x19,0x44,0xc3,0xd2,0x44,0x27,0x41,0x78};
    mont256_from_bytes(&ctx, multiplier, multiplier43);
    mont256_from_bytes(&ctx, multiplicand, multiplicand43);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product43, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier44[32]   = {0x1f,0x26,0x42,0xaa,0xdc,0xde,0xd2,0x04,0x43,0xb3,0x0f,0x66,0x11,0x0e,0x2c,0xb6,0x38,0xef,0xba,0xeb,0xdb,0x31,0xcc,0xd2,0x9b,0xb1,0x83,0xe1,0x15,0x70,0x26,0x6b};
    std::uint8_t multiplicand44[32] = {0xea,0x59,0x67,0x9a,0xed,0x3a,0x32,0xa8,0x6a,0xf2,0x57,0x48,0x8d,0x95,0x9c,0x31,0xfe,0x8a,0xd4,0xa1,0x56,0xd2,0xa6,0x8c,0x02,0xf4,0xb3,0x42,0x74,0x2a,0x80,0x63};
    std::uint8_t product44[32]      = {0x13,0x3b,0x28,0xa3,0xb5,0xd4,0x6b,0x0c,0xa1,0x6b,0x5b,0x7f,0xcd,0x77,0x82,0xcf,0xc7,0xa4,0xf1,0xb6,0x43,0x38,0xbe,0xa7,0x1a,0x96,0x11,0xd9,0xd6,0x05,0x6a,0xd4};
    mont256_from_bytes(&ctx, multiplier, multiplier44);
    mont256_from_bytes(&ctx, multiplicand, multiplicand44);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product44, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier45[32]   = {0xf0,0x29,0x05,0x31,0x3d,0x0a,0x27,0x0b,0xb5,0xa4,0x32,0xcf,0x86,0xe3,0xe7,0x26,0x0b,0x0f,0x87,0x3b,0x21,0x14,0xe0,0x68,0x9f,0x27,0xf5,0x2c,0x44,0x92,0x74,0xd2};
    std::uint8_t multiplicand45[32] = {0xee,0xa7,0xbb,0x64,0x33,0xa7,0x15,0x68,0x2e,0x5f,0x95,0x0c,0x0c,0xe5,0xaf,0x69,0x43,0x0b,0x91,0xed,0x29,0x54,0xba,0x5c,0xf8,0x1e,0x54,0xdd,0x1c,0x05,0x02,0xc6};
    std::uint8_t product45[32]      = {0xb3,0xdc,0xe9,0x06,0x92,0xd6,0xb8,0x85,0xc0,0x59,0x85,0xc4,0xd0,0x26,0xcb,0x3d,0x22,0x3a,0x4b,0xde,0x3a,0x84,0xc8,0xf9,0xd4,0x4f,0x9c,0x13,0x7b,0xb6,0x86,0xd2};
    mont256_from_bytes(&ctx, multiplier, multiplier45);
    mont256_from_bytes(&ctx, multiplicand, multiplicand45);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product45, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier46[32]   = {0x72,0x18,0x88,0xff,0x4a,0x3a,0xdf,0x99,0x34,0xb3,0xff,0x60,0xc2,0x6e,0x7a,0x42,0x87,0xf5,0x3d,0xdd,0x4e,0x14,0xd5,0x71,0xa0,0xf0,0x96,0xda,0x4f,0xde,0xbb,0xec};
    std::uint8_t multiplicand46[32] = {0xfe,0x97,0x7c,0x56,0x04,0xa6,0x56,0x51,0xcd,0xbd,0xe7,0x47,0x58,0xd5,0x0f,0x1b,0x45,0x40,0xf4,0x26,0x2d,0x8a,0xd8,0xc0,0xac,0x12,0x7e,0x93,0x80,0x05,0xce,0x74};
    std::uint8_t product46[32]      = {0x04,0x05,0xcd,0xf5,0xef,0x51,0x68,0x59,0x8a,0x50,0xa3,0x93,0x92,0x13,0xad,0xb4,0x3a,0xe5,0x52,0xf5,0x7a,0x4e,0x77,0x24,0x8c,0x97,0xef,0x16,0x3f,0x86,0xa6,0x8b};
    mont256_from_bytes(&ctx, multiplier, multiplier46);
    mont256_from_bytes(&ctx, multiplicand, multiplicand46);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product46, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier47[32]   = {0xfa,0x61,0x97,0x74,0x8d,0x11,0x8e,0x37,0x81,0x72,0x8a,0x07,0xbb,0xab,0x27,0xf6,0x04,0xb8,0x15,0x7d,0x03,0xed,0xb9,0x20,0x09,0x75,0x83,0x40,0x40,0x1d,0x68,0xfb};
    std::uint8_t multiplicand47[32] = {0xa8,0x87,0xae,0x22,0x1b,0x35,0x41,0x1b,0x72,0x72,0x3b,0x9c,0xef,0x44,0xc0,0xd5,0x3e,0xe4,0xda,0x5a,0x79,0x89,0xe9,0xd0,0x83,0xa4,0xe6,0x29,0x30,0x80,0x38,0x89};
    std::uint8_t product47[32]      = {0xc6,0x11,0x84,0xf4,0x12,0xa8,0x4f,0x39,0x15,0xc1,0x1a,0x13,0xe4,0x56,0xea,0x86,0x1b,0xc0,0x00,0x87,0x6a,0xb4,0xed,0x92,0xa6,0x80,0x6e,0xae,0x3e,0xef,0x05,0x5c};
    mont256_from_bytes(&ctx, multiplier, multiplier47);
    mont256_from_bytes(&ctx, multiplicand, multiplicand47);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product47, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier48[32]   = {0xe3,0x83,0x8b,0x9e,0xd5,0xa9,0x42,0x2a,0x8b,0xc0,0x83,0x11,0x7e,0xb8,0x6c,0x57,0xa8,0x11,0x00,0xa1,0x6e,0xa3,0x30,0xa1,0xa6,0x6d,0x58,0xb5,0xd1,0xa4,0xc0,0x1e};
    std::uint8_t multiplicand48[32] = {0x3a,0xc4,0xda,0x9a,0xfb,0x81,0x39,0x21,0x37,0x16,0x1c,0x16,0xb0,0x0f,0xd7,0xbb,0x4e,0xca,0xde,0xa2,0x81,0xb6,0x2b,0xb5,0xf8,0x66,0x64,0xae,0x64,0xa1,0x49,0xf5};
    std::uint8_t product48[32]      = {0xb5,0xf2,0x94,0x87,0xa7,0x16,0xe8,0x14,0x5f,0x6e,0x83,0x5e,0x56,0xfc,0x78,0x48,0x1c,0xbf,0xf3,0x17,0x28,0xe5,0x43,0x37,0x31,0xd5,0x26,0x09,0x82,0x02,0x27,0x94};
    mont256_from_bytes(&ctx, multiplier, multiplier48);
    mont256_from_bytes(&ctx, multiplicand, multiplicand48);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product48, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier49[32]   = {0x23,0xc4,0x9c,0xae,0xa2,0xcf,0x62,0xba,0xba,0x95,0x88,0x10,0xb4,0xeb,0xf4,0xb6,0xe1,0xc6,0x0a,0xa3,0xd5,0x10,0xbb,0x04,0x32,0xd9,0x0d,0xcd,0x57,0xbb,0x7d,0x97};
    std::uint8_t multiplicand49[32] = {0x03,0xa6,0x39,0x66,0x21,0x3b,0xca,0x7f,0xd6,0x44,0xde,0x2f,0x0d,0xec,0x68,0x23,0xfb,0x5c,0x9d,0x56,0x58,0xf9,0x2d,0xea,0xfd,0x4b,0xd0,0x30,0x67,0x9a,0x44,0xdd};
    std::uint8_t product49[32]      = {0x58,0xd4,0xa2,0xcf,0x47,0x72,0xe9,0x94,0x96,0x54,0x43,0x54,0x98,0xed,0xcb,0x90,0xe1,0xde,0x28,0x88,0xb2,0x0c,0x85,0xb8,0x98,0x44,0xbf,0xac,0x1f,0xc6,0x4a,0x59};
    mont256_from_bytes(&ctx, multiplier, multiplier49);
    mont256_from_bytes(&ctx, multiplicand, multiplicand49);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product49, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier50[32]   = {0x0e,0x2e,0xc4,0x0a,0x29,0xca,0x86,0x2d,0x6e,0x45,0x05,0xf5,0x41,0x6e,0x99,0xb0,0xe1,0x3e,0x21,0x3e,0xbd,0xaa,0xea,0x00,0xa0,0x1d,0x61,0x6f,0x12,0x1a,0xe3,0xe6};
    std::uint8_t multiplicand50[32] = {0xf8,0x8e,0xde,0x10,0xab,0xa8,0xb9,0xb3,0x81,0x85,0x79,0x7c,0xde,0xdb,0x91,0x09,0x61,0x81,0x77,0xff,0xd7,0x5d,0x67,0x69,0xaa,0x4c,0x5c,0x60,0x15,0xa0,0xcc,0xe6};
    std::uint8_t product50[32]      = {0xd3,0xb6,0xc3,0x35,0xd9,0x4b,0xaf,0x9b,0x34,0xa0,0x7f,0x4a,0x7b,0x3c,0xb3,0xe0,0xa7,0x4b,0xb2,0x15,0x4a,0xad,0x94,0x7b,0x0a,0xdf,0x9c,0x34,0xa5,0x61,0xab,0xe5};
    mont256_from_bytes(&ctx, multiplier, multiplier50);
    mont256_from_bytes(&ctx, multiplicand, multiplicand50);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product50, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier51[32]   = {0x2f,0x73,0x3b,0x05,0x75,0x9e,0xb5,0x59,0x0b,0x94,0xaf,0x3a,0x4b,0x05,0xe1,0xae,0xb1,0x53,0xd6,0x9c,0x3e,0x01,0xaa,0xa6,0x99,0x49,0x8a,0xc4,0x48,0x2c,0xc7,0x8e};
    std::uint8_t multiplicand51[32] = {0x54,0x34,0x81,0x56,0xf6,0x37,0xa4,0x68,0x5d,0x38,0x5e,0x06,0x43,0x63,0xe5,0xd9,0x00,0xed,0x6b,0x02,0x72,0x21,0x8f,0xdc,0x44,0xdf,0x96,0xff,0x28,0x54,0x14,0x24};
    std::uint8_t product51[32]      = {0x49,0x34,0x31,0xcd,0xb4,0x29,0x13,0xb7,0x82,0x88,0xfd,0x84,0x68,0x0e,0x9d,0x17,0x5e,0xd5,0xbc,0x09,0x90,0x44,0xfe,0x6a,0x46,0x92,0x33,0xde,0x96,0xc1,0x35,0xb6};
    mont256_from_bytes(&ctx, multiplier, multiplier51);
    mont256_from_bytes(&ctx, multiplicand, multiplicand51);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product51, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier52[32]   = {0xe1,0xe4,0x37,0xb7,0xf7,0x35,0xef,0xe6,0x08,0xd1,0x80,0x11,0x3e,0x94,0x0b,0xb4,0x52,0xd3,0x1e,0x1b,0x8c,0x0d,0x00,0x33,0xfc,0x23,0x25,0xa9,0xf8,0xfd,0xd2,0x08};
    std::uint8_t multiplicand52[32] = {0x15,0x79,0xda,0x0a,0x61,0xb2,0x48,0x0c,0x55,0xd8,0x5e,0x8d,0x00,0x46,0x0d,0x69,0x2e,0xd6,0x54,0x11,0x5b,0x49,0x15,0x61,0x37,0xc6,0x0e,0x98,0x4f,0x3e,0x88,0x5e};
    std::uint8_t product52[32]      = {0xc3,0x1e,0xa6,0xf4,0x94,0x2e,0x7b,0x7e,0xe5,0x4c,0xdc,0xaa,0xa3,0x9b,0x93,0x1c,0x7b,0x7c,0x47,0xa0,0x3d,0x01,0x8b,0xc7,0xca,0x81,0x25,0x96,0xf5,0x0b,0xbe,0xd0};
    mont256_from_bytes(&ctx, multiplier, multiplier52);
    mont256_from_bytes(&ctx, multiplicand, multiplicand52);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product52, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier53[32]   = {0xc6,0xb7,0x89,0xef,0x81,0x36,0x5a,0xcc,0x3f,0x88,0xaf,0x59,0x33,0x73,0x6d,0xcc,0xa7,0xf0,0xc9,0x9e,0x80,0xb5,0x24,0x4a,0x47,0x67,0xe1,0xfa,0x79,0x82,0x3e,0xb2};
    std::uint8_t multiplicand53[32] = {0x96,0x38,0x92,0xa7,0x66,0x46,0x5d,0x28,0x24,0xd4,0x58,0x9c,0x16,0xfa,0x14,0x21,0xd1,0x29,0xd0,0x67,0x43,0xa0,0x8f,0x06,0x17,0x42,0x0e,0x94,0x01,0x44,0x70,0x2b};
    std::uint8_t product53[32]      = {0x74,0x71,0xc7,0x17,0x4e,0x89,0x79,0xf9,0x85,0xa6,0xea,0x07,0xc1,0x88,0x07,0x33,0x30,0x40,0x11,0x4f,0xb9,0x7f,0xeb,0x73,0x97,0x41,0x79,0x2a,0x94,0x73,0x87,0x15};
    mont256_from_bytes(&ctx, multiplier, multiplier53);
    mont256_from_bytes(&ctx, multiplicand, multiplicand53);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product53, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier54[32]   = {0x15,0xa0,0xa8,0xae,0x3b,0x99,0x68,0x70,0xa1,0x32,0x0b,0x9d,0x4d,0xe2,0xf8,0xad,0x4c,0xb5,0x9a,0xa7,0x05,0xc2,0x2d,0x3f,0x64,0xdb,0xc8,0xd3,0x0a,0xaa,0xaf,0x81};
    std::uint8_t multiplicand54[32] = {0xe4,0x8e,0x9e,0x02,0xa8,0x54,0xc8,0x34,0x27,0xbe,0x9a,0xb1,0xc0,0x23,0x6e,0x49,0xda,0x6e,0x6d,0x8e,0x87,0x78,0xf7,0x42,0xf5,0x27,0xb5,0xc2,0x95,0xe8,0xc9,0x3e};
    std::uint8_t product54[32]      = {0x7a,0x80,0x78,0x2e,0x77,0xef,0x7f,0x78,0x7a,0x68,0xde,0x2a,0x59,0xda,0x5a,0x5b,0xbe,0x75,0x05,0x6c,0x20,0x2a,0x50,0xa4,0xa7,0x94,0xc1,0xf4,0xce,0x28,0x6c,0x81};
    mont256_from_bytes(&ctx, multiplier, multiplier54);
    mont256_from_bytes(&ctx, multiplicand, multiplicand54);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product54, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier55[32]   = {0xb8,0x7e,0x4e,0x2b,0x53,0x7d,0x91,0x28,0xc3,0xa9,0xe8,0x89,0x63,0xb7,0x59,0xf5,0x98,0xb8,0x1c,0x66,0xe1,0x0c,0x16,0x7d,0xc8,0xb6,0xea,0xff,0xb7,0x4b,0x58,0x9b};
    std::uint8_t multiplicand55[32] = {0x25,0x0e,0x7b,0x34,0xa4,0xaa,0x07,0xb4,0x9e,0x63,0x97,0xd4,0xb9,0x62,0x45,0xd3,0x48,0xbf,0xcb,0xcf,0x26,0x43,0x37,0x98,0x7e,0x83,0x49,0x04,0xfc,0x17,0x34,0x98};
    std::uint8_t product55[32]      = {0xa8,0xd2,0x63,0x83,0x73,0x7e,0x9b,0x39,0x5f,0xb1,0x5f,0x78,0x97,0x98,0xda,0xe1,0x06,0x52,0xd9,0x96,0x65,0xe9,0xe8,0x43,0xfc,0xfa,0x98,0x68,0x27,0xa6,0x56,0xa1};
    mont256_from_bytes(&ctx, multiplier, multiplier55);
    mont256_from_bytes(&ctx, multiplicand, multiplicand55);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product55, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier56[32]   = {0x6d,0xe2,0xfb,0x1f,0xa0,0x98,0xd6,0x91,0x83,0x52,0xbc,0x85,0xe4,0x56,0x55,0x9c,0xb7,0x0a,0xf5,0xf2,0xd5,0xd5,0x89,0x1f,0xd3,0x29,0xd6,0x5c,0x0b,0x35,0xb1,0xde};
    std::uint8_t multiplicand56[32] = {0xc0,0xbb,0xe6,0xed,0x86,0x14,0xf5,0x04,0xe8,0xee,0x65,0xa1,0x23,0xa9,0xa9,0xda,0x81,0x6b,0x23,0x32,0xcf,0xed,0x94,0x3b,0xb3,0x78,0x3a,0x7c,0xbb,0xdd,0xbb,0x9b};
    std::uint8_t product56[32]      = {0x7e,0xee,0x5f,0x7f,0x1e,0x33,0x58,0x1d,0x68,0xe4,0x18,0xa0,0x9c,0x00,0xc4,0x8e,0xaf,0x2c,0x5a,0xee,0xdc,0x29,0xc3,0x17,0x09,0x38,0x4b,0xc0,0x3f,0x19,0x0f,0x01};
    mont256_from_bytes(&ctx, multiplier, multiplier56);
    mont256_from_bytes(&ctx, multiplicand, multiplicand56);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product56, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier57[32]   = {0xaf,0xbc,0x9c,0xa9,0xd3,0x8f,0x8c,0x45,0x04,0x1d,0xcd,0x94,0xcd,0xff,0x5a,0x1c,0xd0,0x1a,0x91,0x4c,0xd5,0xbe,0x78,0x5a,0x91,0x87,0xdf,0x42,0x81,0x1e,0x76,0x16};
    std::uint8_t multiplicand57[32] = {0xa4,0x94,0x6d,0x15,0xb1,0x7d,0xd2,0x55,0xf4,0xc1,0x82,0x26,0xae,0xd2,0x3b,0x0f,0xb6,0x10,0x4b,0x84,0xe4,0x90,0x7d,0x49,0xcc,0x47,0x93,0xd7,0x95,0x85,0x0e,0x21};
    std::uint8_t product57[32]      = {0x7e,0x25,0x82,0xc4,0xee,0x43,0x6a,0xa6,0xac,0x86,0xff,0x91,0x05,0x44,0xa8,0x8a,0xf8,0xa6,0xd7,0x4d,0x3e,0x48,0x83,0x4d,0xcf,0x04,0x81,0x57,0x9a,0xa1,0x99,0x42};
    mont256_from_bytes(&ctx, multiplier, multiplier57);
    mont256_from_bytes(&ctx, multiplicand, multiplicand57);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product57, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier58[32]   = {0xf5,0xa2,0xd8,0x79,0x5c,0x57,0x53,0x2b,0xa3,0x1a,0x49,0xdd,0x22,0x12,0x65,0x40,0x0a,0xb7,0x79,0x88,0x07,0xfa,0x22,0xf7,0x15,0xc8,0x91,0xff,0x3a,0xdd,0x65,0x27};
    std::uint8_t multiplicand58[32] = {0x04,0xd2,0xbe,0x09,0xa0,0xb5,0x58,0x64,0x0c,0xff,0xf0,0x54,0x8e,0xfb,0xa4,0x42,0x73,0x8e,0x0b,0x77,0xd5,0xf8,0x60,0xc3,0x60,0x6a,0x0d,0xeb,0x1a,0xdb,0xce,0x5d};
    std::uint8_t product58[32]      = {0xf9,0x92,0x04,0x96,0x28,0x7d,0x51,0xc0,0x77,0x3c,0xf2,0x4b,0x77,0xda,0x72,0x17,0x68,0x31,0xa1,0x4a,0xaf,0x01,0x3f,0xf7,0xab,0x75,0x42,0x98,0xcd,0xe3,0x36,0xdd};
    mont256_from_bytes(&ctx, multiplier, multiplier58);
    mont256_from_bytes(&ctx, multiplicand, multiplicand58);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product58, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier59[32]   = {0x74,0xfa,0x94,0x12,0x00,0xd9,0x35,0x34,0x43,0x87,0xee,0x7b,0x7d,0x42,0x64,0x6f,0x3e,0x9b,0x76,0x8f,0xae,0x40,0x01,0xe3,0x88,0x0c,0xb4,0x01,0xa0,0x50,0x60,0x98};
    std::uint8_t multiplicand59[32] = {0x17,0x89,0x81,0x9f,0x89,0x02,0xda,0xfc,0xe5,0xd9,0xfe,0x81,0x80,0xc2,0xb5,0xf1,0xee,0xb8,0x9f,0xf1,0xbf,0x8e,0x51,0xaa,0x11,0xf2,0xd4,0x4d,0xcc,0x35,0xe8,0x34};
    std::uint8_t product59[32]      = {0x69,0x23,0x82,0xc3,0x64,0x3d,0x84,0xff,0x14,0x07,0xe9,0x0c,0x37,0x92,0x11,0x2c,0x4d,0x23,0xbe,0x9b,0x48,0xb3,0xad,0x0f,0x8f,0x64,0x55,0x71,0xec,0x27,0xc5,0x5f};
    mont256_from_bytes(&ctx, multiplier, multiplier59);
    mont256_from_bytes(&ctx, multiplicand, multiplicand59);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product59, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier60[32]   = {0xcf,0x28,0xf6,0x5e,0x40,0x8f,0xc1,0x46,0x79,0x4e,0xc9,0x26,0xbc,0x9e,0x28,0xea,0xbe,0xe8,0x06,0x26,0x10,0xe8,0xad,0x01,0x86,0xa7,0x4a,0x63,0xa8,0xc7,0xd9,0xe0};
    std::uint8_t multiplicand60[32] = {0x3b,0x11,0x85,0xd9,0x34,0x89,0x22,0xd7,0xc1,0xa6,0x24,0xdc,0xba,0xb5,0xb3,0x73,0x3c,0x1a,0xe9,0x17,0x43,0xfb,0x9f,0xbc,0xd8,0x9c,0x36,0xb2,0x13,0x0f,0x27,0xb2};
    std::uint8_t product60[32]      = {0x1e,0x36,0x8f,0x41,0x3c,0x55,0x34,0xad,0xed,0x32,0x6a,0x81,0xbb,0xf3,0xff,0x85,0xb2,0xb9,0x58,0x0f,0xe2,0x35,0x16,0x5f,0x6d,0x25,0xdc,0xdb,0x47,0x5d,0xf4,0x1b};
    mont256_from_bytes(&ctx, multiplier, multiplier60);
    mont256_from_bytes(&ctx, multiplicand, multiplicand60);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product60, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier61[32]   = {0x13,0xa5,0x39,0x7f,0x61,0xef,0x7b,0xd1,0xd8,0x74,0xbc,0x79,0x7e,0x73,0x6d,0x5f,0x75,0xd8,0xd8,0xa4,0xf9,0xc9,0xc6,0x79,0xa6,0x61,0xf6,0x2c,0xbd,0x65,0x68,0x0c};
    std::uint8_t multiplicand61[32] = {0xa1,0xfe,0xb6,0x24,0x9d,0xf2,0x02,0x5f,0x0b,0xf7,0xa4,0xbd,0xc4,0x58,0x27,0x2f,0x49,0x8d,0xbf,0xa8,0xaf,0x06,0xbc,0xf7,0xe9,0x14,0x57,0xdb,0x7a,0xa0,0x68,0xf1};
    std::uint8_t product61[32]      = {0xe2,0x35,0x4b,0x71,0xfe,0xbf,0xce,0x58,0x41,0x0a,0x3d,0x22,0x98,0x9b,0x89,0x24,0x51,0xb7,0x58,0x1d,0x98,0xd9,0x2f,0x34,0x00,0x0e,0x2a,0x04,0x32,0xe6,0xf9,0xd0};
    mont256_from_bytes(&ctx, multiplier, multiplier61);
    mont256_from_bytes(&ctx, multiplicand, multiplicand61);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product61, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier62[32]   = {0xa6,0xca,0xf4,0xa3,0x41,0x02,0x3a,0xed,0x54,0xef,0x12,0x5a,0x25,0xbd,0xa6,0x59,0x99,0x86,0x48,0xe0,0x13,0xd5,0x31,0x6f,0x32,0xc3,0x24,0x44,0xa4,0x8c,0x1d,0x5c};
    std::uint8_t multiplicand62[32] = {0x7b,0x7f,0xec,0x4b,0x03,0x31,0x2e,0xad,0x22,0x29,0x30,0xae,0x91,0x58,0xd4,0xa8,0x9f,0x03,0xbc,0x5a,0x4d,0xee,0x48,0x12,0xb1,0x61,0x07,0xf1,0xbe,0x43,0x7c,0x7b};
    std::uint8_t product62[32]      = {0xbc,0xc1,0x0b,0x57,0x0a,0x24,0xe2,0x88,0xd6,0xf7,0x8c,0x14,0xec,0xc2,0x58,0x23,0xee,0x17,0x0d,0x48,0xd9,0x3b,0xf8,0x68,0x67,0xe8,0x4d,0x69,0x67,0xe4,0xbc,0xa9};
    mont256_from_bytes(&ctx, multiplier, multiplier62);
    mont256_from_bytes(&ctx, multiplicand, multiplicand62);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product62, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier63[32]   = {0x37,0xba,0xc2,0x33,0xb1,0x33,0x0c,0x3f,0x19,0x7a,0x14,0xe2,0xac,0x08,0x4b,0xa5,0xf8,0xf6,0x59,0xac,0x44,0xce,0x4a,0xb3,0x7c,0x5d,0x42,0xdc,0x0f,0x87,0x7a,0xe3};
    std::uint8_t multiplicand63[32] = {0x77,0x45,0x10,0xca,0x76,0xf4,0x25,0x1e,0x49,0x19,0x61,0xa1,0x84,0x3b,0xae,0xe9,0xb5,0x78,0x90,0x9c,0x4a,0x75,0x91,0xf2,0x7d,0x57,0x5d,0x17,0xac,0xfb,0x2d,0x5e};
    std::uint8_t product63[32]      = {0x61,0xfc,0x7e,0x11,0xae,0xdd,0xfd,0xba,0x1f,0x4e,0x70,0xb3,0xd3,0x00,0x05,0x8f,0x88,0x32,0x7d,0x60,0x04,0x63,0x01,0x94,0x31,0x9d,0x86,0x8f,0x0c,0x31,0xe6,0x87};
    mont256_from_bytes(&ctx, multiplier, multiplier63);
    mont256_from_bytes(&ctx, multiplicand, multiplicand63);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product63, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier64[32]   = {0x4f,0xc9,0xe9,0x18,0x33,0x02,0x0c,0xcd,0x8c,0x90,0x47,0x3e,0xe4,0xc7,0x17,0xfd,0xfe,0x48,0xef,0x63,0x1e,0x56,0x34,0x08,0xc4,0x65,0x3c,0xde,0x77,0x62,0x00,0xb5};
    std::uint8_t multiplicand64[32] = {0x13,0x93,0x29,0x04,0x75,0x7f,0x1c,0xba,0x4a,0x22,0x7f,0x39,0x04,0x7b,0x2c,0x10,0x79,0x12,0xef,0x4a,0xef,0xae,0x5d,0x4e,0x15,0xfa,0x8b,0x65,0xfa,0x66,0x72,0xcd};
    std::uint8_t product64[32]      = {0x76,0x76,0xf6,0xab,0xc5,0x9a,0x0c,0x58,0x12,0xd1,0x06,0x03,0x6a,0xc4,0x09,0xf3,0xd7,0x01,0x79,0xfc,0x32,0x2f,0x1e,0xec,0xe7,0x0d,0x57,0x3a,0x54,0xde,0x09,0xc4};
    mont256_from_bytes(&ctx, multiplier, multiplier64);
    mont256_from_bytes(&ctx, multiplicand, multiplicand64);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product64, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier65[32]   = {0x63,0x08,0x7e,0x52,0x44,0xc6,0xb8,0x95,0xfe,0x74,0x9e,0x67,0x73,0x0f,0x37,0xf1,0xfe,0x9e,0xb4,0xad,0xf7,0xd5,0xf1,0x24,0x81,0xb1,0xc0,0x25,0xd1,0xe4,0xd0,0xa3};
    std::uint8_t multiplicand65[32] = {0x17,0x1e,0x1a,0x8c,0x94,0xdb,0x5f,0x8f,0x13,0x19,0xd4,0x24,0x35,0xf1,0x03,0x00,0xee,0x37,0x9c,0x65,0xf2,0x12,0x01,0xe4,0xea,0xa3,0x55,0x6c,0x35,0xb7,0xe4,0x48};
    std::uint8_t product65[32]      = {0xcd,0x7b,0x08,0x47,0x82,0xc9,0x0c,0x03,0x7f,0xe6,0x31,0xf7,0xcd,0x73,0x35,0x17,0x4c,0x75,0x69,0xf2,0x8f,0x62,0xec,0x94,0x3c,0x7c,0xb1,0xc4,0x46,0xca,0x9f,0x51};
    mont256_from_bytes(&ctx, multiplier, multiplier65);
    mont256_from_bytes(&ctx, multiplicand, multiplicand65);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product65, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier66[32]   = {0x9a,0x76,0x2d,0x54,0x21,0xf2,0x67,0xe2,0x5c,0x0b,0xb4,0x0f,0xf3,0xe6,0xca,0x73,0x43,0x05,0xe9,0x86,0x86,0x29,0x2b,0xb5,0xbf,0x5b,0x41,0x1b,0x24,0x49,0x1d,0xf6};
    std::uint8_t multiplicand66[32] = {0x5d,0x7c,0xfe,0xd1,0xb4,0x0d,0xe5,0x6d,0x1c,0xd8,0x6f,0xc1,0xe3,0x09,0x66,0x19,0x47,0x91,0xc2,0xe9,0x82,0x3d,0x11,0xed,0xa1,0xb5,0x01,0xd6,0xd1,0xf9,0xbd,0xfe};
    std::uint8_t product66[32]      = {0x44,0x1b,0x35,0xab,0x4a,0xf8,0xe3,0x81,0xa4,0x72,0x2c,0x73,0x7c,0xc9,0x15,0xff,0xd8,0x6a,0xdf,0x59,0x8c,0x18,0xaa,0x0b,0x4e,0xd6,0x83,0x05,0x35,0x1f,0x6d,0xaf};
    mont256_from_bytes(&ctx, multiplier, multiplier66);
    mont256_from_bytes(&ctx, multiplicand, multiplicand66);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product66, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier67[32]   = {0x28,0xb8,0x80,0x73,0x06,0x5b,0x8c,0x35,0x64,0xe2,0x76,0x02,0x7c,0x73,0xb6,0xc9,0xe0,0x4b,0x0d,0xce,0xe5,0xd0,0x0a,0x4d,0x7f,0x75,0x95,0xb5,0x3b,0x3b,0xf4,0xbf};
    std::uint8_t multiplicand67[32] = {0xba,0x28,0xa6,0x79,0x4d,0x4c,0xa9,0xc7,0x67,0xc9,0x8f,0xb9,0x73,0x65,0x06,0xec,0xae,0x7c,0x8f,0x09,0x7d,0xdf,0xcb,0xc9,0xf3,0x30,0x8c,0xe5,0x00,0xeb,0x4e,0x11};
    std::uint8_t product67[32]      = {0xed,0xfb,0x79,0x10,0xb3,0xd8,0xf4,0x66,0xf3,0x03,0x63,0xdc,0x91,0xa4,0x4a,0xbe,0x3f,0xea,0xa6,0x05,0x95,0xf2,0x1b,0x64,0xf9,0x35,0x8a,0xc2,0x03,0x45,0x0c,0xdc};
    mont256_from_bytes(&ctx, multiplier, multiplier67);
    mont256_from_bytes(&ctx, multiplicand, multiplicand67);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product67, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier68[32]   = {0x54,0xd1,0xac,0x6b,0xd7,0x19,0x61,0x89,0x1e,0xf3,0xea,0x44,0x50,0xea,0x7d,0xa7,0x60,0x48,0x7e,0x15,0x58,0x0d,0xc5,0xab,0x6a,0x8a,0xd9,0xcb,0x24,0x05,0x63,0x60};
    std::uint8_t multiplicand68[32] = {0xf0,0x9c,0x0a,0xfb,0x1e,0xbb,0x07,0x94,0x65,0xf4,0x56,0xaa,0xd6,0xcf,0xf7,0x18,0x56,0x99,0x08,0xf6,0xc0,0x30,0x1b,0x21,0x53,0x15,0x8c,0xe4,0x00,0x72,0x1f,0x84};
    std::uint8_t product68[32]      = {0x0d,0x5d,0x10,0x60,0xd3,0xb6,0x85,0x2e,0x0b,0xf8,0x45,0x45,0xce,0x58,0x67,0x78,0x23,0x29,0xf1,0x59,0xc5,0xd5,0x68,0xa8,0xdb,0x0d,0xb2,0xff,0xe0,0xff,0x72,0xc1};
    mont256_from_bytes(&ctx, multiplier, multiplier68);
    mont256_from_bytes(&ctx, multiplicand, multiplicand68);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product68, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier69[32]   = {0x40,0xd2,0x84,0x06,0x4a,0x32,0x7e,0x2d,0xbd,0x6a,0x99,0x6d,0xe6,0xcd,0x10,0xf1,0x03,0x00,0x30,0x05,0xb6,0x88,0xb6,0x61,0x32,0x1c,0x17,0x44,0xed,0x28,0x79,0xc1};
    std::uint8_t multiplicand69[32] = {0x13,0x8e,0xfe,0xf9,0x96,0xd4,0x48,0x0f,0xde,0xb6,0x7a,0xe7,0xff,0xb0,0xdd,0x9e,0x63,0xe1,0x98,0x69,0x64,0x95,0x0d,0xc2,0x10,0xa2,0x5b,0x19,0x5f,0x49,0xf0,0xfc};
    std::uint8_t product69[32]      = {0xdd,0x56,0x93,0x5f,0x92,0xe0,0x35,0xf6,0xd7,0xb5,0xfd,0xff,0xb6,0xac,0x37,0xe2,0x40,0xee,0x31,0xb8,0xc3,0xb9,0x4e,0x4a,0x0b,0xa0,0x06,0xba,0xa0,0x4d,0x7c,0x81};
    mont256_from_bytes(&ctx, multiplier, multiplier69);
    mont256_from_bytes(&ctx, multiplicand, multiplicand69);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product69, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier70[32]   = {0x47,0xd7,0xdf,0x79,0x0c,0x5b,0x4c,0x59,0xda,0xb0,0x79,0x29,0x46,0x70,0x93,0x12,0xc1,0x72,0xb2,0x98,0x6d,0x94,0xdd,0x6d,0xec,0xe8,0x07,0x99,0x5c,0x57,0x72,0x2e};
    std::uint8_t multiplicand70[32] = {0x26,0x1f,0x40,0xdf,0xef,0x82,0xd1,0xa3,0xa2,0x8c,0xf7,0xb1,0x49,0x1e,0x99,0xf5,0xa9,0x77,0x66,0xfb,0xd5,0xad,0x53,0x60,0x0d,0x36,0xce,0x2c,0x1a,0x09,0xa8,0x40};
    std::uint8_t product70[32]      = {0x33,0xc5,0x64,0x90,0x62,0x33,0x0f,0x26,0x53,0xe2,0x1d,0x85,0x70,0xbe,0xde,0x2c,0xb6,0x1f,0x54,0xd5,0x22,0x16,0xf8,0x0d,0xef,0x60,0x14,0xd0,0x72,0x1e,0x1f,0x05};
    mont256_from_bytes(&ctx, multiplier, multiplier70);
    mont256_from_bytes(&ctx, multiplicand, multiplicand70);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product70, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier71[32]   = {0xc5,0xef,0x5c,0xfb,0x30,0x99,0xf2,0x71,0x50,0xcb,0x40,0x7a,0x82,0xce,0x78,0x6f,0x6f,0xad,0x79,0x36,0x44,0x06,0xc0,0x53,0xf8,0x95,0xfc,0x55,0x3f,0xd3,0xbe,0x98};
    std::uint8_t multiplicand71[32] = {0xc2,0xfb,0xd8,0xa3,0xcf,0xdc,0xc2,0x57,0x07,0x6d,0x49,0x0a,0xe2,0x5f,0x4b,0x1c,0x6d,0x80,0xde,0x7c,0xf4,0xc7,0x3f,0x2b,0xc8,0xff,0x1c,0x38,0x5f,0x93,0xd1,0x80};
    std::uint8_t product71[32]      = {0xa7,0xb4,0xb8,0xe6,0x2b,0x1d,0x19,0x48,0xbd,0x84,0x0e,0x6e,0x31,0x27,0x8e,0xa8,0x50,0x63,0x27,0xde,0x11,0x7a,0xc2,0xc1,0xac,0x91,0x0d,0xf8,0x67,0xf6,0x71,0x2d};
    mont256_from_bytes(&ctx, multiplier, multiplier71);
    mont256_from_bytes(&ctx, multiplicand, multiplicand71);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product71, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier72[32]   = {0x34,0x14,0x5e,0x87,0x8c,0x9a,0x37,0x51,0x8d,0xdc,0xf8,0x3c,0xf0,0xd1,0xab,0x56,0xe0,0x2f,0x9a,0x72,0xe9,0xd6,0x25,0xc9,0x66,0x69,0x21,0x58,0xa1,0x82,0x63,0x27};
    std::uint8_t multiplicand72[32] = {0x9d,0x6b,0x02,0x3f,0x73,0x6b,0x96,0xa0,0x69,0x2f,0xd3,0x60,0xbb,0x7b,0x73,0x8e,0xee,0xf7,0x95,0xcd,0x0c,0xaa,0x76,0x12,0x14,0xa0,0xb0,0x0b,0xb8,0x35,0xe8,0xa5};
    std::uint8_t product72[32]      = {0x61,0xae,0x69,0x3f,0xe5,0x2b,0x7f,0xb5,0x15,0xfd,0x04,0x42,0xc9,0x89,0x6a,0x94,0x75,0xf9,0x8b,0xde,0x65,0xd8,0xd8,0x57,0x7f,0xfa,0xc5,0x46,0x43,0xd4,0x02,0xd4};
    mont256_from_bytes(&ctx, multiplier, multiplier72);
    mont256_from_bytes(&ctx, multiplicand, multiplicand72);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product72, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier73[32]   = {0xe9,0x72,0x9f,0x3f,0x0c,0x89,0xc0,0x01,0x7c,0x4e,0xa6,0x03,0x49,0x44,0xf2,0xce,0xde,0x96,0x2a,0x6d,0xa4,0xfd,0x57,0xc5,0x23,0x79,0x7d,0x45,0xc0,0xae,0xd9,0xc5};
    std::uint8_t multiplicand73[32] = {0x48,0x20,0x82,0x31,0x57,0xfa,0x49,0xe5,0x6a,0x34,0xb3,0x71,0x78,0xe1,0x0e,0x70,0x2b,0xb7,0x1c,0x68,0x20,0x97,0x79,0x8c,0x8c,0xd3,0xe4,0x18,0xed,0x41,0x42,0xba};
    std::uint8_t product73[32]      = {0xcd,0xe5,0x8c,0xd5,0xd7,0x83,0x86,0xe6,0x57,0x26,0xda,0xd9,0xf6,0x40,0x06,0x13,0xe1,0x76,0xb0,0xf0,0xbf,0x86,0x63,0x81,0xf8,0xaa,0x58,0x2c,0xe0,0x4b,0x8b,0x3b};
    mont256_from_bytes(&ctx, multiplier, multiplier73);
    mont256_from_bytes(&ctx, multiplicand, multiplicand73);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product73, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier74[32]   = {0x67,0xfd,0x54,0x99,0x42,0x9a,0x70,0x79,0xa7,0x1f,0x11,0xb2,0xf9,0xee,0x8b,0xc8,0xbd,0x1e,0x69,0x12,0xbd,0x31,0x3b,0xee,0x41,0x78,0x5b,0xc6,0x4c,0x3a,0xc6,0xfc};
    std::uint8_t multiplicand74[32] = {0x1e,0xa7,0x72,0x28,0x64,0xf5,0x49,0x69,0xab,0x3b,0x74,0xfe,0x8e,0xac,0xa2,0x88,0x7b,0xb1,0xd1,0x24,0x4d,0x03,0x9b,0x72,0x3d,0x19,0x26,0xac,0xa7,0xef,0x4f,0x5d};
    std::uint8_t product74[32]      = {0x02,0xc4,0x33,0xb5,0xea,0x10,0x18,0x88,0xdc,0xa5,0x5f,0x3d,0xbe,0x2c,0xf9,0xda,0x92,0x1f,0xcc,0x9e,0xde,0xc4,0x1c,0x38,0x41,0xab,0x95,0x77,0x80,0xe7,0x16,0x06};
    mont256_from_bytes(&ctx, multiplier, multiplier74);
    mont256_from_bytes(&ctx, multiplicand, multiplicand74);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product74, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier75[32]   = {0xcf,0xd3,0xdd,0x72,0xe7,0xec,0xfd,0x0c,0x80,0x27,0xa2,0xa2,0x35,0x37,0x22,0x35,0x13,0x3e,0x61,0x53,0x29,0x62,0x59,0xc8,0xa4,0xa9,0x15,0xd0,0x2a,0xd6,0x4c,0xe9};
    std::uint8_t multiplicand75[32] = {0xc2,0x5e,0x11,0x4f,0xff,0x18,0xfe,0x33,0x55,0x34,0xa0,0x34,0xe8,0x00,0x9d,0x90,0x73,0xf6,0xe5,0x3d,0x38,0x53,0x93,0x3d,0x8c,0xe6,0x21,0xef,0x7f,0x40,0x5b,0xc8};
    std::uint8_t product75[32]      = {0xdb,0xf2,0xff,0xc2,0x05,0x9d,0xb9,0x57,0xfe,0x21,0x03,0x28,0xda,0x00,0x8a,0xca,0xdc,0x5a,0x37,0x2e,0x88,0xbb,0x88,0x1a,0x74,0x49,0x33,0xfb,0xf0,0x3e,0x32,0xc4};
    mont256_from_bytes(&ctx, multiplier, multiplier75);
    mont256_from_bytes(&ctx, multiplicand, multiplicand75);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product75, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier76[32]   = {0x2c,0xb8,0xd1,0x4c,0x17,0x39,0x10,0xe3,0x3e,0x7c,0x65,0x67,0x31,0x41,0x97,0x75,0x8c,0x3b,0xa8,0x59,0x23,0xbc,0x91,0x52,0x6d,0x6b,0x98,0x7a,0x73,0x30,0x9b,0x95};
    std::uint8_t multiplicand76[32] = {0xcf,0x32,0x1d,0x63,0x42,0x23,0xb8,0xaa,0x5e,0x49,0x42,0x2a,0x3d,0x37,0x66,0x42,0x51,0xbc,0xd7,0x7a,0x17,0x51,0xf5,0x79,0x8e,0x4d,0xc3,0xa3,0x57,0x8a,0x60,0xd8};
    std::uint8_t product76[32]      = {0x15,0x89,0xce,0x4c,0xc5,0xb4,0xed,0x1d,0x9d,0x62,0xd2,0x0e,0xea,0x1c,0xf5,0x88,0x87,0xa1,0x51,0x0f,0xbf,0x0e,0xfd,0x4f,0x2b,0xa0,0x0b,0xd6,0x81,0x32,0x40,0xad};
    mont256_from_bytes(&ctx, multiplier, multiplier76);
    mont256_from_bytes(&ctx, multiplicand, multiplicand76);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product76, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier77[32]   = {0x62,0x01,0xa9,0xd3,0x69,0xac,0x0f,0x03,0xde,0xe0,0xa8,0x43,0xbf,0xe9,0x8f,0x8c,0x05,0x24,0x13,0x7f,0xe3,0x22,0xe9,0x6d,0x33,0xbf,0x91,0x57,0x91,0xd2,0x77,0xf2};
    std::uint8_t multiplicand77[32] = {0xc0,0x8a,0x58,0xd7,0x56,0x94,0x7a,0x7a,0x45,0x2e,0x70,0x4d,0x60,0x7a,0x47,0x32,0x35,0xc2,0xe2,0x29,0x86,0x2f,0xe2,0x31,0xbe,0xef,0x67,0xfb,0x69,0xf4,0x46,0x12};
    std::uint8_t product77[32]      = {0x34,0xc8,0x27,0x07,0x34,0x64,0x54,0xc1,0x64,0x6f,0xa5,0x36,0xf6,0x8c,0x8e,0xe8,0xa4,0xf1,0x5d,0xc9,0x1b,0xd0,0x46,0x49,0x59,0x8c,0x79,0x08,0xe5,0x94,0x8e,0xee};
    mont256_from_bytes(&ctx, multiplier, multiplier77);
    mont256_from_bytes(&ctx, multiplicand, multiplicand77);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product77, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier78[32]   = {0xaf,0xcf,0x0e,0x77,0x20,0x39,0x43,0xf6,0x5c,0x32,0x7a,0x6d,0xf7,0xba,0x38,0xb6,0x93,0x04,0x10,0x6e,0x47,0x0b,0x4f,0xad,0x7f,0x86,0x7d,0x5f,0x0f,0xe3,0x21,0xec};
    std::uint8_t multiplicand78[32] = {0x17,0xb4,0x83,0x4c,0x37,0x49,0x5c,0x5e,0xd9,0x3f,0xf7,0x16,0xdc,0xe4,0x7b,0x21,0xca,0x51,0xe1,0x52,0xa1,0x2f,0x3a,0x94,0x87,0x7b,0x55,0xcb,0x80,0xde,0x8b,0x3e};
    std::uint8_t product78[32]      = {0x15,0x26,0xa5,0x6f,0x44,0x63,0xc6,0x3c,0x1c,0xb5,0x3f,0xcb,0xef,0xe0,0xe7,0x7c,0x44,0x9d,0x18,0xc2,0x94,0x77,0x95,0x64,0xe3,0x6e,0x8a,0x6a,0xb9,0x94,0xaa,0x24};
    mont256_from_bytes(&ctx, multiplier, multiplier78);
    mont256_from_bytes(&ctx, multiplicand, multiplicand78);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product78, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier79[32]   = {0x6e,0x8c,0xd9,0x4e,0x72,0x23,0xc6,0x8a,0xa5,0x52,0x9b,0x05,0x66,0x56,0x7b,0xc4,0x62,0x72,0x92,0xf8,0x3f,0x9a,0xa8,0x84,0xe5,0x94,0x09,0xc1,0x45,0x61,0x9f,0xc0};
    std::uint8_t multiplicand79[32] = {0x20,0x93,0x42,0xca,0x05,0x95,0x5f,0xb9,0xf7,0xd1,0x7e,0xbd,0xdf,0x75,0xc8,0x83,0xd0,0x78,0x84,0xb7,0xd9,0x43,0x55,0x41,0x4f,0xe0,0x48,0x02,0xf4,0x35,0xa5,0x73};
    std::uint8_t product79[32]      = {0xe5,0xc4,0x1b,0x0f,0xbc,0xf1,0xb3,0xa1,0xc6,0xd1,0xaf,0x32,0x65,0xaf,0xe7,0x92,0xdb,0x2d,0x32,0xb8,0x10,0x41,0xc0,0x78,0x35,0x90,0x28,0x54,0x73,0xb3,0x9a,0x22};
    mont256_from_bytes(&ctx, multiplier, multiplier79);
    mont256_from_bytes(&ctx, multiplicand, multiplicand79);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product79, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier80[32]   = {0xf7,0xe1,0x47,0xfd,0x79,0x28,0x1c,0x19,0xcd,0xe3,0x47,0xab,0xe5,0x4c,0x5d,0xe6,0xc3,0x81,0x3c,0xe6,0xb5,0xa2,0x90,0x61,0x6c,0xd9,0xe6,0x2a,0x08,0x41,0x1c,0x07};
    std::uint8_t multiplicand80[32] = {0xed,0x9b,0xf0,0xb6,0xed,0x44,0x8d,0x4e,0xee,0x24,0x1c,0x43,0x64,0x3a,0xb9,0xe2,0x12,0xb9,0x2a,0x01,0x00,0x0b,0xb5,0xf9,0x7d,0x65,0x21,0x35,0x96,0x51,0x32,0xd6};
    std::uint8_t product80[32]      = {0x16,0xc9,0xdf,0x3f,0xe7,0x0a,0x45,0xc2,0xbb,0x8f,0x56,0x3a,0x4c,0x13,0xab,0xc7,0x54,0x01,0xdb,0x31,0x8f,0x8a,0xc5,0x7b,0xec,0x30,0xe4,0x77,0x15,0x7e,0x91,0xdf};
    mont256_from_bytes(&ctx, multiplier, multiplier80);
    mont256_from_bytes(&ctx, multiplicand, multiplicand80);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product80, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier81[32]   = {0xc8,0x79,0xb6,0x63,0x3f,0x9b,0x6b,0xb2,0x72,0xee,0x6a,0x2e,0xf8,0xe4,0xcb,0x5c,0x77,0xd8,0xc5,0x69,0xda,0xff,0x9a,0x0b,0x87,0x21,0xec,0xf8,0xd3,0x59,0xd0,0x7a};
    std::uint8_t multiplicand81[32] = {0x1b,0xe0,0x3d,0xf0,0xae,0x9c,0x78,0xbd,0xf8,0xcd,0x9e,0xc3,0x85,0xb9,0xc0,0x9a,0x26,0xed,0xf1,0xbd,0x27,0x85,0x57,0x98,0x39,0x4a,0xfb,0xe9,0x1b,0xea,0x70,0x5e};
    std::uint8_t product81[32]      = {0xbe,0x43,0x9e,0xa5,0x3c,0x99,0xcd,0xb8,0x21,0x17,0xc4,0x30,0x4c,0x81,0x8f,0x7d,0xb6,0x26,0x6e,0x56,0x43,0xc5,0x85,0xff,0x23,0x44,0x5e,0xda,0x0b,0x3f,0x72,0x08};
    mont256_from_bytes(&ctx, multiplier, multiplier81);
    mont256_from_bytes(&ctx, multiplicand, multiplicand81);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product81, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier82[32]   = {0xe5,0x17,0x4e,0xbd,0xc3,0xc9,0xf7,0xe3,0xd8,0xb4,0xc8,0x31,0xa5,0xb8,0x9b,0x2f,0xb3,0x74,0xfa,0xb6,0xb8,0xc3,0xa4,0xd2,0xd3,0x4d,0x1c,0x0d,0xf1,0x05,0x86,0x67};
    std::uint8_t multiplicand82[32] = {0x20,0x2a,0xb6,0xfa,0xc8,0x44,0xb8,0xfd,0x00,0x59,0x86,0x5a,0x0a,0x1f,0xb4,0x3b,0xc6,0xe0,0x67,0x3a,0x8d,0x2f,0x29,0xe7,0x15,0xc2,0xc8,0x1a,0x75,0x13,0x41,0x07};
    std::uint8_t product82[32]      = {0xd4,0x6f,0x40,0x2e,0x6f,0x26,0x43,0x78,0xb8,0x66,0x2b,0x07,0x3f,0xc7,0x76,0xd0,0xa6,0x19,0x6d,0x06,0x5d,0xa0,0x21,0x3f,0xee,0xa7,0x68,0x96,0x3c,0x73,0x69,0xdc};
    mont256_from_bytes(&ctx, multiplier, multiplier82);
    mont256_from_bytes(&ctx, multiplicand, multiplicand82);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product82, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier83[32]   = {0xf6,0x62,0x22,0x2e,0x4d,0xc4,0xac,0x8c,0xb7,0x0b,0xa8,0x58,0xa5,0x3f,0xdd,0xc9,0x09,0x9f,0x9c,0x9f,0xeb,0x7f,0xe2,0x6b,0x91,0xc3,0x09,0x8c,0x3b,0x8a,0x27,0xba};
    std::uint8_t multiplicand83[32] = {0xc3,0x8b,0x48,0xa2,0xb2,0xd6,0x43,0xa2,0x6f,0xfb,0x72,0x6a,0xa2,0xe3,0xf9,0x3a,0x87,0x3b,0x99,0x03,0x40,0x75,0x91,0x6e,0xa0,0x60,0x84,0x6c,0x20,0xc2,0x6f,0x71};
    std::uint8_t product83[32]      = {0x6f,0xfe,0xef,0x60,0x3e,0x5b,0x82,0x3e,0x62,0x57,0x9d,0x49,0xfd,0xf8,0xa8,0x80,0x2a,0x88,0xcd,0x78,0xcb,0xb4,0x6c,0x47,0x54,0x74,0xa6,0xf0,0xe3,0x15,0x73,0x80};
    mont256_from_bytes(&ctx, multiplier, multiplier83);
    mont256_from_bytes(&ctx, multiplicand, multiplicand83);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product83, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier84[32]   = {0x31,0x13,0x5d,0xe9,0x95,0x38,0x57,0xd7,0xf1,0x8b,0xde,0x0e,0x86,0x41,0x7b,0x60,0x4c,0xe3,0xb0,0xcc,0x12,0x02,0x95,0x2f,0x19,0x75,0x36,0xb1,0x1c,0xb4,0xba,0x55};
    std::uint8_t multiplicand84[32] = {0x89,0x98,0x0c,0x50,0x02,0xad,0x9d,0x2b,0x00,0x4b,0x7f,0xd0,0x99,0xdf,0x20,0x9b,0xca,0x5d,0x5e,0x7d,0x39,0x3c,0xbc,0xdd,0x42,0xc9,0x27,0xb9,0x63,0x59,0x56,0xbe};
    std::uint8_t product84[32]      = {0xac,0x02,0xfd,0x65,0xef,0x52,0x53,0x60,0x04,0xf0,0x82,0xf0,0xb7,0x08,0x4b,0xa0,0x5a,0xe5,0xfc,0x54,0xaa,0xad,0xa5,0xb6,0x76,0x62,0x56,0x0a,0x3a,0xe5,0x8c,0x0b};
    mont256_from_bytes(&ctx, multiplier, multiplier84);
    mont256_from_bytes(&ctx, multiplicand, multiplicand84);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product84, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier85[32]   = {0xd6,0xe3,0xa7,0x1e,0xa5,0x02,0xe8,0xa8,0x50,0xfc,0xc6,0x26,0xf5,0x7d,0x17,0x09,0x47,0x52,0x91,0x94,0x75,0xef,0xd2,0x33,0xff,0x12,0x5e,0xb4,0x4d,0x30,0x7f,0xe4};
    std::uint8_t multiplicand85[32] = {0x07,0x7e,0xf3,0x2a,0x3f,0x3f,0x37,0xea,0x8c,0x08,0x56,0xa4,0x3c,0x19,0xc3,0x15,0x86,0xba,0x22,0xdd,0x79,0xad,0x89,0x99,0x3e,0x0b,0x25,0xcd,0xe2,0x3f,0x03,0xcc};
    std::uint8_t product85[32]      = {0x83,0xae,0x3b,0x22,0x0e,0x38,0x3b,0x80,0x50,0x1e,0xab,0x16,0x77,0x64,0x8f,0xe4,0x1d,0x69,0x3e,0x63,0xb0,0x88,0x2d,0xe8,0x67,0x13,0x55,0x2d,0xfa,0x57,0x35,0xe1};
    mont256_from_bytes(&ctx, multiplier, multiplier85);
    mont256_from_bytes(&ctx, multiplicand, multiplicand85);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product85, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier86[32]   = {0x31,0xb1,0x89,0x1a,0x05,0x93,0xdb,0xa2,0x0e,0x28,0xb6,0x4f,0x4e,0xb1,0x9f,0xca,0xa6,0x4f,0x76,0x13,0xb4,0x64,0x2e,0xa4,0x69,0x6c,0x63,0xd6,0xf5,0xea,0xd0,0x65};
    std::uint8_t multiplicand86[32] = {0x3a,0x53,0xc1,0x76,0x41,0xdb,0x89,0x8e,0x14,0xc2,0x73,0x2a,0x6b,0x86,0x29,0x0b,0xa5,0xac,0xd3,0x41,0xac,0xa9,0x9f,0xd0,0xe2,0x85,0x6e,0xc6,0x7f,0x91,0x42,0x86};
    std::uint8_t product86[32]      = {0xb5,0x5c,0x8d,0x97,0xa1,0x77,0xd9,0xfa,0x02,0x1c,0xd4,0xde,0x52,0x1c,0xe3,0xb0,0x67,0x51,0xad,0x40,0x4f,0xd3,0xa3,0xb8,0x8d,0xe0,0x82,0x5b,0x71,0x7b,0x60,0x79};
    mont256_from_bytes(&ctx, multiplier, multiplier86);
    mont256_from_bytes(&ctx, multiplicand, multiplicand86);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product86, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier87[32]   = {0xb2,0x21,0x71,0x39,0x08,0xba,0x9b,0xd9,0x7e,0x31,0x8a,0xd6,0x3a,0x0e,0xa6,0xe1,0x5e,0xc6,0x9b,0xe3,0xec,0xd7,0x57,0x0b,0x6c,0xa0,0x64,0x96,0xaa,0xd7,0xc7,0xc0};
    std::uint8_t multiplicand87[32] = {0x01,0xba,0x98,0x5a,0x32,0xb5,0x58,0xfd,0x65,0x77,0xbb,0x54,0xae,0xbc,0xb0,0xaa,0x5c,0xc0,0xff,0x06,0x6b,0xa9,0x9d,0x01,0xb7,0xe4,0x9f,0x36,0x56,0x8a,0x8c,0x29};
    std::uint8_t product87[32]      = {0xfb,0x68,0x8d,0x2c,0xe6,0xf3,0x2f,0xf0,0x0c,0xc6,0x14,0x95,0xfa,0xb2,0x5e,0x77,0xc9,0x67,0x79,0xd4,0x65,0xd0,0x53,0xd2,0xd3,0xa0,0x66,0xdc,0xee,0xa4,0xda,0x4f};
    mont256_from_bytes(&ctx, multiplier, multiplier87);
    mont256_from_bytes(&ctx, multiplicand, multiplicand87);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product87, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier88[32]   = {0x7e,0xe5,0xe8,0x57,0x34,0x89,0x34,0x98,0x11,0x43,0x40,0xff,0x81,0x3f,0xb5,0xcd,0xd8,0x5b,0xbb,0x6b,0xbd,0x37,0x92,0x9d,0x4a,0xc7,0xcc,0xc3,0xcc,0x0c,0x66,0x82};
    std::uint8_t multiplicand88[32] = {0x77,0x11,0xb7,0x57,0x3b,0x16,0x49,0x43,0x31,0xa5,0x9c,0x4a,0xd1,0xeb,0xd0,0x86,0xc4,0x0f,0x36,0x09,0x4f,0xcc,0x9a,0x5c,0x33,0x4e,0x51,0xaf,0xf8,0x48,0xa9,0x56};
    std::uint8_t product88[32]      = {0x62,0x47,0x86,0x3a,0x81,0x49,0xd7,0x28,0x0c,0xa9,0x53,0xce,0x44,0x58,0xa6,0xed,0x68,0xec,0x38,0x21,0x03,0x26,0x87,0x2f,0xd8,0xe5,0x4d,0xd2,0x79,0x3d,0xa0,0xc5};
    mont256_from_bytes(&ctx, multiplier, multiplier88);
    mont256_from_bytes(&ctx, multiplicand, multiplicand88);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product88, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier89[32]   = {0x9f,0xa4,0x0d,0xd6,0xf3,0xb1,0x7a,0xf0,0x1b,0xe7,0xf3,0xcf,0x4b,0x80,0xb8,0x28,0xe3,0xab,0x62,0x83,0xc2,0xae,0x35,0xd2,0x43,0xd8,0x7a,0x97,0x38,0xb0,0x79,0xe1};
    std::uint8_t multiplicand89[32] = {0xe9,0x0f,0xb6,0x51,0x6a,0xc2,0x6a,0xe0,0x7c,0x2c,0x6a,0x87,0x39,0x2b,0xc5,0x52,0xe5,0x7f,0x76,0x91,0x2f,0xf3,0xc2,0x3c,0x9c,0x2f,0x67,0x23,0x7e,0xea,0x6f,0xe1};
    std::uint8_t product89[32]      = {0xa5,0x1d,0xbb,0x42,0xcd,0xa8,0x0b,0xbc,0x92,0xb7,0x31,0x9b,0x85,0x72,0xa5,0xa1,0x14,0x15,0x5e,0x9b,0x8d,0x5a,0x71,0x9f,0x9c,0xfa,0xab,0x6d,0x0d,0x08,0xe7,0xb9};
    mont256_from_bytes(&ctx, multiplier, multiplier89);
    mont256_from_bytes(&ctx, multiplicand, multiplicand89);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product89, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier90[32]   = {0x0d,0xea,0x6e,0x4e,0x64,0xb9,0xcb,0x1c,0xec,0x03,0x2e,0x6b,0x25,0x79,0x5c,0x18,0x98,0x44,0xf4,0x76,0xf2,0xe2,0x05,0x4d,0x0e,0x71,0x59,0x7a,0xaa,0x50,0xb9,0x6f};
    std::uint8_t multiplicand90[32] = {0xb5,0xb9,0x4a,0xf3,0x0d,0x45,0x6b,0xe0,0x6a,0x56,0xaa,0xc3,0x24,0x54,0x48,0xc8,0x98,0x9b,0xc9,0xdc,0xf9,0x5f,0xe8,0xa0,0x06,0x0c,0x88,0x04,0x36,0x83,0xd4,0xbc};
    std::uint8_t product90[32]      = {0x4d,0x7e,0x72,0x4f,0x63,0x0a,0xf3,0xb6,0x39,0xdc,0xab,0x86,0x7d,0x90,0xaf,0xf2,0x2c,0x59,0xff,0x18,0xcd,0x6d,0xd1,0x0f,0x37,0x9c,0x43,0x16,0x01,0xfc,0x3e,0x98};
    mont256_from_bytes(&ctx, multiplier, multiplier90);
    mont256_from_bytes(&ctx, multiplicand, multiplicand90);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product90, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier91[32]   = {0x50,0x6f,0x68,0xac,0xe2,0x32,0x89,0x94,0xb6,0x47,0xe8,0xa8,0xe5,0xee,0x4c,0x91,0x73,0x1b,0xbc,0x41,0x64,0xb0,0xbb,0x14,0x2f,0x21,0x7e,0x72,0x0f,0x65,0x06,0x38};
    std::uint8_t multiplicand91[32] = {0x30,0xd0,0xa2,0xb8,0x54,0x49,0x40,0xe1,0x2a,0x66,0xf9,0x13,0xee,0x7d,0x0a,0xe2,0x14,0x51,0x03,0xc7,0xff,0x5e,0x1d,0x1f,0x1c,0xfb,0x0a,0x06,0xbb,0x93,0xc8,0xeb};
    std::uint8_t product91[32]      = {0x9b,0x57,0x98,0x4d,0x91,0x5d,0x74,0x98,0x43,0x3d,0xc9,0xbd,0x9a,0x2a,0xc2,0x27,0x46,0x88,0x92,0x63,0x12,0x3a,0x56,0x07,0x57,0x3f,0xbd,0x7e,0x84,0x6f,0x88,0x97};
    mont256_from_bytes(&ctx, multiplier, multiplier91);
    mont256_from_bytes(&ctx, multiplicand, multiplicand91);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product91, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier92[32]   = {0x4f,0xd3,0xe7,0x58,0x08,0x2a,0x2f,0x4d,0x77,0xb5,0xab,0xcb,0xbf,0x0e,0x11,0xe0,0x86,0x59,0x22,0x43,0xef,0x95,0xee,0xe8,0xa7,0x08,0x28,0xa7,0x2f,0x7d,0xba,0x08};
    std::uint8_t multiplicand92[32] = {0x71,0x43,0x6e,0x1d,0x54,0xea,0x20,0x61,0xfc,0x27,0xd6,0x83,0x5f,0xb6,0xd6,0x25,0xd6,0xd1,0x06,0xfb,0x60,0xed,0x33,0xa0,0xb9,0xb2,0x53,0xe3,0xaa,0x18,0x13,0x45};
    std::uint8_t product92[32]      = {0xc2,0x85,0xd3,0xca,0x39,0xdc,0xfa,0x4b,0x3b,0xda,0xf3,0x81,0xd0,0x83,0x1e,0x64,0x31,0x0d,0xb3,0xc2,0xac,0x56,0x61,0x7e,0x91,0x00,0x65,0xf9,0x95,0xaa,0xff,0x73};
    mont256_from_bytes(&ctx, multiplier, multiplier92);
    mont256_from_bytes(&ctx, multiplicand, multiplicand92);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product92, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier93[32]   = {0x6b,0x91,0x1f,0x97,0x59,0xf9,0xbb,0x79,0x14,0xac,0xe1,0xcb,0x47,0xa1,0x64,0xe4,0x14,0x07,0xab,0x33,0x00,0xbc,0x22,0xcb,0x1b,0xe4,0xa5,0xdb,0x2b,0x54,0xaf,0x77};
    std::uint8_t multiplicand93[32] = {0x61,0x50,0x2d,0xee,0x35,0x18,0x53,0x76,0xc2,0x41,0x0a,0xd1,0xf6,0xda,0x7a,0x63,0x8f,0xa6,0x24,0xf7,0x1f,0xab,0x58,0x84,0xe2,0x9a,0xac,0xea,0xf4,0x9c,0x9e,0xba};
    std::uint8_t product93[32]      = {0xba,0xbf,0x55,0xe2,0x9f,0x53,0x86,0xb8,0x04,0x3c,0x76,0xa0,0x32,0xb8,0xa6,0x71,0x2a,0xa5,0x66,0x7a,0x0c,0x7b,0xd6,0x57,0x91,0xd1,0x46,0xb5,0xc7,0x48,0xcd,0x8c};
    mont256_from_bytes(&ctx, multiplier, multiplier93);
    mont256_from_bytes(&ctx, multiplicand, multiplicand93);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product93, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier94[32]   = {0x16,0x77,0x74,0xef,0x6e,0xb4,0xff,0xf8,0xcd,0xce,0xc4,0x08,0xd2,0x6f,0x1d,0x76,0x4f,0x06,0xe9,0x5a,0xd2,0x52,0xa6,0x17,0xc4,0xcb,0xa0,0x38,0x5b,0x4c,0x0d,0x73};
    std::uint8_t multiplicand94[32] = {0x72,0x43,0xd4,0x7c,0xeb,0x64,0xc5,0xc4,0x8a,0xa1,0xa5,0x9c,0x5f,0x6a,0x35,0xd9,0x32,0x1a,0x6e,0xc1,0x79,0x34,0xf0,0xb8,0xb4,0x8b,0xb0,0x75,0x0c,0x9c,0x20,0xef};
    std::uint8_t product94[32]      = {0xa4,0xe8,0x2e,0xe2,0x79,0x05,0xd5,0xdf,0xd4,0x2f,0x48,0xb5,0x13,0xec,0x04,0x56,0x98,0x08,0x79,0x34,0x22,0x06,0xa4,0xb5,0xdd,0xe4,0xf0,0x55,0x5c,0x07,0x52,0x74};
    mont256_from_bytes(&ctx, multiplier, multiplier94);
    mont256_from_bytes(&ctx, multiplicand, multiplicand94);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product94, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier95[32]   = {0xa1,0xb4,0x9b,0xf7,0x07,0xc0,0x90,0x9c,0x79,0x7b,0x15,0x38,0xe5,0xa1,0x5b,0x79,0xbc,0xc0,0xfd,0x98,0x5d,0x3f,0x69,0xce,0x52,0xc4,0x64,0x1b,0x31,0x6a,0x2a,0x12};
    std::uint8_t multiplicand95[32] = {0x60,0x25,0x33,0xdc,0x0a,0x68,0x01,0x3d,0x67,0x9f,0x2d,0x9e,0xc4,0x44,0x5a,0xae,0xa0,0x1a,0xc2,0x3a,0xcf,0xd3,0xbb,0x74,0x3f,0x7d,0xc8,0x6b,0x69,0x2a,0x4f,0x0e};
    std::uint8_t product95[32]      = {0x6b,0x35,0x63,0x3c,0x07,0xea,0x73,0x4f,0x82,0x87,0x6f,0x24,0x1a,0xb2,0x8c,0xfd,0xe6,0x05,0x72,0x28,0x03,0xb7,0x26,0x43,0x56,0xeb,0x04,0x44,0x2f,0x69,0xec,0xd9};
    mont256_from_bytes(&ctx, multiplier, multiplier95);
    mont256_from_bytes(&ctx, multiplicand, multiplicand95);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product95, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier96[32]   = {0x31,0xe7,0xae,0xd1,0x41,0xcb,0xcc,0x3a,0x0f,0xdf,0x7c,0xc6,0xeb,0x8a,0x25,0xfc,0xcd,0xa7,0x90,0x77,0x10,0x05,0x3d,0x2c,0x76,0xcc,0x05,0x73,0x08,0xec,0x37,0x9a};
    std::uint8_t multiplicand96[32] = {0x55,0xc0,0xa7,0x4d,0x45,0xb6,0x69,0xf7,0x5c,0xeb,0xe2,0x13,0x56,0xcd,0x42,0xd2,0x9b,0x09,0xab,0x55,0xe6,0x07,0x7d,0x79,0x10,0x17,0x0d,0x2b,0xbf,0x4e,0x30,0x2c};
    std::uint8_t product96[32]      = {0x5d,0x21,0xa1,0x79,0x0b,0x9b,0xd4,0x8f,0xe5,0x9b,0x6d,0x51,0xe6,0xc2,0x32,0x1d,0xef,0x2b,0x54,0xe1,0x18,0x9e,0xb3,0x0f,0xd5,0x97,0xc0,0x12,0x1e,0x43,0xbc,0xb9};
    mont256_from_bytes(&ctx, multiplier, multiplier96);
    mont256_from_bytes(&ctx, multiplicand, multiplicand96);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product96, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier97[32]   = {0xb0,0x88,0x24,0x11,0xb7,0x75,0x70,0xa4,0xbf,0x16,0x8d,0xa7,0x43,0x1d,0xbc,0x3f,0x0b,0x28,0x6c,0x70,0x9d,0xf2,0x4d,0x5e,0xf4,0x29,0xc6,0x22,0xf5,0x2b,0x25,0x49};
    std::uint8_t multiplicand97[32] = {0x98,0x77,0x27,0x90,0xc1,0x72,0x6f,0x06,0xb8,0xb8,0xf2,0x70,0x00,0xf7,0x2d,0x3c,0x4c,0x22,0xca,0xb7,0x46,0x8f,0xb5,0x96,0xec,0x9a,0x36,0x0c,0x51,0x05,0x12,0x2a};
    std::uint8_t product97[32]      = {0x4a,0x55,0xec,0xdd,0x56,0x99,0xe9,0xda,0x31,0xc3,0xbb,0xf9,0x65,0x46,0x15,0xf9,0x27,0x34,0x74,0x05,0xcb,0xfc,0x81,0xc6,0x41,0xd5,0x6a,0x67,0xfa,0xd6,0x0e,0x2a};
    mont256_from_bytes(&ctx, multiplier, multiplier97);
    mont256_from_bytes(&ctx, multiplicand, multiplicand97);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product97, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier98[32]   = {0xd3,0x75,0xef,0xf1,0x06,0x35,0xaf,0xef,0x10,0xb9,0x9a,0xc9,0xf1,0x78,0xd7,0x7f,0xf2,0x4d,0x04,0xfd,0xa2,0x4c,0x84,0x07,0xce,0x3f,0xa0,0x28,0xea,0x9d,0x18,0xb2};
    std::uint8_t multiplicand98[32] = {0xc6,0xbf,0x4f,0xa2,0xf4,0x33,0x7b,0xd1,0x77,0x3a,0xfe,0x02,0xf4,0xef,0x61,0x42,0xb7,0x2f,0xac,0x4a,0x79,0xa5,0xfd,0x62,0x1b,0x75,0x7b,0x20,0x3b,0xde,0xa8,0xc3};
    std::uint8_t product98[32]      = {0xe2,0x5f,0x62,0xe9,0x42,0x04,0x57,0xee,0x8d,0x65,0x62,0x31,0x1e,0x6a,0x82,0xca,0xcd,0xec,0x98,0xa6,0x1d,0x77,0x0d,0x47,0xdf,0x1a,0xa7,0x39,0x48,0x56,0x18,0xcc};
    mont256_from_bytes(&ctx, multiplier, multiplier98);
    mont256_from_bytes(&ctx, multiplicand, multiplicand98);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product98, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
    std::uint8_t multiplier99[32]   = {0x21,0xf9,0x1a,0x99,0x7e,0x54,0x4d,0x56,0xd0,0x96,0xbf,0xd6,0x6e,0x10,0x6c,0x0e,0xe9,0xde,0x04,0x79,0x40,0x44,0x9a,0xa0,0xca,0x30,0x42,0x18,0x62,0xf2,0xa2,0x1b};
    std::uint8_t multiplicand99[32] = {0x4d,0xa6,0x09,0x90,0xbd,0x0d,0x8c,0xfe,0xee,0x59,0xb3,0x97,0xcd,0x75,0x1e,0x08,0x02,0x3a,0x80,0xa2,0x2e,0xd5,0x1b,0x12,0x7f,0x1d,0x49,0x0e,0xed,0x97,0xec,0x76};
    std::uint8_t product99[32]      = {0x2f,0x58,0xb9,0xad,0xeb,0xa4,0x2e,0xf6,0x50,0xc5,0x97,0xd8,0x90,0xf8,0xdc,0x05,0xe6,0x8f,0x03,0xe6,0xa5,0xa8,0x1c,0xcd,0x07,0x91,0x12,0x02,0x5f,0x57,0xa0,0x11};
    mont256_from_bytes(&ctx, multiplier, multiplier99);
    mont256_from_bytes(&ctx, multiplicand, multiplicand99);
    mont256_mul(&ctx, product, multiplier, multiplicand);
    mont256_to_bytes(&ctx, product_data, product);
    if (std::memcmp(product99, product_data, 32) != 0)
    {
        throw std::runtime_error("err in mont256_mul");
    }
}
